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Abstract 



A real spin-charge separation scheme is found based on a saddle-point state 



fSl ■ of the t — J model. In one-dimensional (ID) case, such a saddle-point repro- 

o 

duces the correct asymptotic correlations at the strong-coupling fixed-point of 
the model. In two-dimensional (2D) case, the transverse gauge field confining 

a 

i . spinon and holon is shown to be gapped at finite doping so that a spin-charge 

deconfinement is obtained for its first time in 2D. The gap in the gauge fluc- 

O ■ 

tuation disappears at half- filling limit, where a long-range antiferromagnetic 
^ . order is recovered at zero temperature and spinons become confined. The 

most interesting features of spin dynamics and transport are exhibited at fi- 
nite doping where exotic residual couplings between spin and charge degrees 
of freedom lead to systematic anomalies with regard to a Fermi-liquid sys- 
tem. In spin dynamics, a commensurate antiferromagnetic fluctuation with 
a small, doping-dependent energy scale is found, which is characterized in 
momentum space by a Gaussian peak at (n/a, tt /a) with a doping-dependent 
width (oc y/S, 5 is the doping concentration). This commensurate magnetic 
fluctuation contributes a non-Korringa behavior for NMR spin-lattice relax- 
ation rate. There also exits a characteristic temperature scale below which a 
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pseudo-gap behavior appears in the spin dynamics. Furthermore, an incom- 
mensurate magnetic fluctuation is also obtained at a finite energy regime. In 
transport, a strong short-range phase-interference leads to an effective holon 
Lagrangian which can give rise to a series of interesting phenomena including 
linear-T resistivity and T 2 Hall-angle. We discuss the striking similarities of 
these theoretical features with those found in the high-T c cuprates and give a 
consistent picture for the latter. Electronic properties like Fermi-surface and 
superconducting pairing in this framework are also discussed. 
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I. INTRODUCTION 



The normal state of high-T c cuprate superconductors has shown peculiar properties 
in both charge and spin channels. In transport aspect, the resistivity^ exhibits a linear- 
temperature dependence up to 1000 K and down to a temperature ~ T c , which can be 
as low as 10 K. This temperature dependence has been related to a scattering rate! that 
behaves like ri^-T with rj ~ 2 for all the optimally doped cuprates with T c ranging from 
10 up to over 100 K (a linear-frequency dependence of the scattering rate at frequency 
> ksT/h is also found in the infrared spectroscopy^ up to 0.15 eV). Such a linear-T longi- 
tudinal resistivity is also accompanied by a Hall coefficient! which implies hole-like charge 
carriers and shows a 1/T dependence in contrast to a T- independent Fermi liquid case. 
A Hall-angle experimental analysis!!! has demonstrated that the 1/T behavior in the Hall 
coefficient is due to an additional scattering rate in the transverse channel, which behaves 
like T 2 . Most recently, the magneto-resistance has been found to have a T -4 temperature 
dependence,! at variance with the Kohler's rule. Furthermore, the thermopower has shown 
a strong doping-dependence,0 which decreases with the increase of doping and even changes 
sign in the overdoped regime. All of these transport results are anomalous with regard to 
the canonical phenomena in a conventional Fermi liquid (FL) system. 

Spin magnetic properies have also exhibited a number of anomalies, which persist into 
the superconducting phases. Two powerful probes of spin dynamics in the cuprates are 
nuclear-magnetic-resonance! (NMR) and neutron-scattering.! In the cuprates, NMR spin- 
lattice relaxation rate Tf -1 , which probes a very small energy scale (~ 10~ 4 meV), has shown 
a strong enhancement! at low temperature for planar Cu nuclei, in contrast to the usual 
Korringa behavior in a Fermi liquid. On the other hand, T-f 1 for planar O nuclei! is more 
or less close to the Korringa law. The sharp contrast of the NMR T{~ between the planar 
Cu and O nuclcei strongly suggest sliiHl! that the non-Korringa behavior at Cu sites should 
be caused by antiferromagnetic (AF) correlations among the Cu spins whose effect could be 
effectively cancelled out at O sites. AF spin fluctuations have been also directly measured 
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by the inelastic neutron-scattering experiment si at higher energies in La%~ x [Sr, Ba] x CuO^ 
(LSCO) [Ref. and YBa 2 Cu 3 7 - S (YBCO) [Ref. materials. However, in con- 

trast with more or less universal behaviors in NMR spin relaxation rates, neutron-scattering 
data have shown distinctive characteristics among these compounds. For example, commen- 
surate AF fluctuations are observed in the insulating LSCO systems00 and underdoped 
YBCO compounds^lS with characteristic energy scales much smaller than the exchange 
energy J ~ 120 meV. But incommensurate AF correlations have been founds in the metal- 
lic LSCO down to an energy scale as low as 1 meV. Nevertheless, recent analyses have 
revealediiii that the NMR relaxation rates in LSCO are inconsistent with an incommen- 
surate AF structure at u ~ 0. This is because the effect of incommensurate fluctuations 
could not be well cancelled out at O sites, and a large non-Korringa signal would leak to 
the latter, as in contradiction to experiment. In the optimally-doped YBCO, the absence of 
low-energy AF fluctuation in neutron-scattering experiment!^ is also inconsistent with the 
NMR measurement s.i Thus a consistent explanation of the NMR data and a reconciliation 
of NMR and neutron-scattering are two essential issues in spin dynamics. 

Given these anomalous transport and magnetic properties, one would naturally question 
if a single type of elementary excitations carrying both charge and spin can explain them 
without leading to any intrinsic contradiction. Since the sharp Fermi edge has been well 
identified in the cuprates by the angle-resolved photo-emission measurements,@@ such low- 
lying electron-like excitations, if exist, should always locate near the Fermi surface, no 
matter whether it is a Fermi liquid or a marginal Fermi liquid (MFL) [Ref. system. 
The main difference between FL and MFL lies in their energy spaces (the broadening of 
the quasi-particle in MFL is linear in energyi! as compared to u 2 in FL case). But in the 
phase space (momentum space), they are facing essentially the same problems.^ Low-lying 
states as labelled by momenta are uniformly distributed in momentum space around Fermi 
surface. The only conceivable structure would come from the topology of Fermi surface. 
The question is whether such a structure of Fermi surface (maybe with nesting and van 
Hove singularity) is capable of explaining an army of abundant experimental anomalies. In 



transport channel, one of the key problems is how to get hole-type charge carriers, since the 
elementary excitation is of electron-type. Even though a special curvature of Fermi surface 
may explain the hole-like sign indicated by the Hall effect, the total charge carrier number 
could by no means coincide with the hole number as measured with regard to the half-filled 
insulating parent compound. Many other problems could arise within such a framework. 
For example, a monotonic decrease of the thermopoweiQ from over lOOfjY/K down to and 
even changing sign, when doping is increased from zero to an overdoped level, is difficult to 
understand here because the phase-space is too limited. Even with the presence of a van Hove 
singularity, an overall change as estimated by theory (~ 8fiV/K [Ref. [28| ) is too small as 
compared to the experimental value. Without additional structure, a second scattering rate 
impliedSi in the Hall-angle measurements is even more difficult to comprehend, especially 
the different ways by which it gets into the transverse and longitudinal channels. 

In the spin channel, a number of FL and MFL theoriesil@ have been proposed in terms 
of the geometry of Fermi surface. The topology of Fermi surface determines essential char- 
acteristics of spin dynamics through the Lindhard-type response function. These theories 
provide an explanation of incommensurate AF structure^ in the metallic LSCO, but fail to 
account for other equally important issues. One is about the NMR relaxation rates whose 
aforementioned distinctive behaviors at planar Cu and O nuclei are difficult to reconcile .HH 
The fundamental reason is that a Lindhard response function cannot produce a sharp enough 
AF peak near (ir/a, n/a). Furthermore, the overall width of the AF peak in these theories 
is generally not sensitive to doping concentration, and one cannot find NMR anomalies at 
low temperature to be sufficiently enhanced with the decrease of doping as indicated by the 
experiments.! Most seriously, the itinerant picture fails to provide a small AF energy scale, 
as the only characteristic scale is the Fermi energy ep. Although a peak at small energy 
may be introduced by van Hove singularity]^ it is hard to be related to a magnetic energy. 
The exchange energy scale ~ J as indicated in various magnetic measurements, and an even 
smaller doping-dependent magnetic energy scale found by inelastic neutron-scatteringBH^ 
is a strong indication that the local spin description of AF fluctuation may be more ap- 



propriate than an itinerant picture. This is further supported by the high temperature Tf" 1 
measurement in LSCO [Ref. |3(| which shows that doping effect is basically evaporated when 
T > where T^ 1 at finite doping coincides with that of an insulating antiferromagnet, 

described by the localized spins under Heisenberg model.!]] 

Therefore, the phase-space limitation in an itinerant description, which becomes in- 
evitable in a system where fermionic elementary excitations carry both spin and charge, 
prevents a consistent explanation of rich phenomena in the high-T c cuprates. This implies 
that an electron in the cuprates may be a composite particle, consisting of more elementary 
excitations which carry spin and charge quanta (spinon and holon) separately. In this way, 
spinon and holon may find enough phase-spaces of their own to account for experimental 
anomalies. A spin-charge separation scenario for the cuprates has been first proposed by 
Anderson. HH Under this idea, a strong on-site Coloumb repulsion can lead to a split of 
energy band with a large Hubbard gap, and thus real charge carriers (holons) in the con- 
duction band will have equal number as that of doped holes controlled by experiments. It is 
a significant first step towards understanding the transport properties. However, in order to 
reconcile with a sharp Fermi edge found by the angle-resolved photo-emission, a fermionic 
spinon with a similar Fermi surface has to be assumed. @ Thus at least in the spin channel, 
basic properties should not be very different from conventional FL or MFL theory, which is 
a serious set back. Later it has been further realizediHil that there actually exists a strong 
coupling among spinons and holons, known as gauge interaction, and such a gauge force 
plays a role essentially to confined spinon and honon together, like the quark confinementi! 
in quantum chromodynamics. In other words, whether there is a real spin-charge separation 
in a 2D strongly-correlated model is still unclear. 

So we are in a paradox situation. Experimental measurements, on one hand, point to 
a possibility of spin-charge separation. Strong-correlation theories, on the other hand, lead 
to spin-charge confinement. This awkward situation is actually caused by an oversimplified 
choice to let spinons to have a large Fermi surface, which undermines a local description 
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for spins. One may argue that in order to recover a correct electron Fermi surface-an 
experimental constraint-a large spinon Fermi surface would be necessary. But this is not the 
case even in one dimension, where a real spin-charge separation exists in the Hubbard model 
and an electron Fermi surface (points) still satisfies^ the Luttinger volume theorem. Even 
though usually people may be used to the thinking^ that a spinon Fermi surface is the reason 
for an electron Fermi surface in ID, the various correlations have never been correctly derived 
under such a picture. In fact, a FL description of spinons is a too rough approximation in 
ID. In the strong-coupling (large U) regime, a correct spin-charge separation descriptionil 
has been established in a path-integral formalism,^ where an electron is described as a 
composite particle of a spinon and a holon, together with a nonlocal phase-shift field. It is 
this phase-shift field that helps to recover the right Fermi surface position. In this formalism, 
spinons are described by a local spin representation without a Fermi surface, and both spin 
and density, as well as various electronic correlation functions have been correctly obtained.^ 
An important lesson learned here is that the electronic Fermi surface is no longer essential 
in a spin-charge separation scheme, and it can be reproduced from a phase-shift in a correct 
spin and charge separation scheme. (We note that in ID an another powerful method is 
the bosonization method&"0 in weak-coupling regime. There the non-Fermi-liquid like low- 
lying spin and charge processes can be directly described near the electronic Fermi points. 
But its possible generalization to two dimensions has many difficulties and is still under 
investigation.) 

This leads to a new 2D spin-charge scenario for the cuprates. In this scenario, spinons 
would be free of "duty" to be responsible for electron Fermi surface, so that they could 
get sufficient phase-space freedom to describe anomalous spin dynamics. Electronic Fermi 
surface would be produced by extra phase-shift field associated with the fermion-statistics of 
electron, and should not be directly related to holon and spinon. This provides a fundamental 
reason for the Luttinger theorem to be valid even in strong correlations. It also means a 
Fermi-surface topology is no longer crucial in determining spin and charge dynamics. That 
is, the shape of a Fermi surface due to the detailed band structure, which may vary from 



one material to others in the cuprates, should not be so relevant to the basic spin and charge 
anomalies in a spin-charge separation framework. Therefore it may well justify using a simple 
one-band t — J model@'@^ to describe the realistic cuprates if a spin-charge separation is 
indeed present there. It is noted that several interesting experimental features have been 
already known for the t — J model. For example, spin degree of freedom is described in a 
local spin representation, which reduces to the Heisenberg Hamiltonian at half-filling and 
well describes0ij] the magnetic behaviors in the insulating cuprates. On the other hand, 
charge carriers in this system are naturally holes as measured from the half-filled insulating 
phase, which is also a strong indication of the experimental relevance for the model. A 
spin-charge separation is already well-known for such a model in ID case. 

It is the purpose of the present paper to establish the above-described spin-charge sep- 
aration scheme within the t — J framework. Like other approximations, we cannot directly 
prove that this 2D spin-charge separation state is the solution of the t — J model. But such 
a state will satisfy the following important criteria. The transverse gauge fluctuation will be 
found suppressed in the long-wavelength and low-energy regime so that spinon and holon 
are indeed deconfined at finite doping (i.e., spin-charge separation). This is in contrast to 
singular infrared gauge fluctuationil~Il in the case of uniform RVB state. At zero doping 
limit, the gap of the gauge fluctuation will vanish such that the spinons become confined 
again to form spin-wave excitations. In this case, a long-range AF order will emerge at zero 
temperature. Furthermore, the present state can also reproduce the correct results in ID 
case, which is an important check of the theory because only in ID one has exact solution.!! 

In the present paper, the effective Hamiltonian describing the spin-charge separation 
saddle-point will be derived in the following form: 



Heff = H s + H h , 




where H s and are the spinon and holon Hamiltonians, respectively, defined by 



H. 



s 





<ij><y 



s 



H h = -t h ^ (S~*% +a U) hfh 3 + h.c. (1.3) 

<ij> 

Here bi a and hi are called spinon and holon annihilation operators, respectively, which are 
connected to electron operator in Eq.(1.6) below. H s and Hh in Eqs.(1.2) and (1.3) would 
resemble a standard tight-binding Hamiltonian if there are no phase fields A^-, and 
■A-ij (< ij > refers to two nearest-neighboring sites on the lattice). Here [as defined 
in Eq.(2.45)] represents a 7r-flux threading through each plaquette. A tight-binding model 
under can be easily diagonalized.il The most prominent feature in Eqs.(1.2) and (1.3) is 
the presence of A? 7 - and A^-, the so-called topological phases. A^ is defined by 

(1.4) 

with complex coordinate Zi — Xi + iyi in 2D, and nf = hfhi as the holon number operator. 
Ajj is given by 

A-,- = - V Im In 

with n\ a = b^bia- as the spinon number operator. 

The usual gauge coupling, which otherwise would confine spinon and holon together, 
will be shown suppressed and has been neglected in Eq.(l.l). Thus the residual interaction 
between spinon and holon are solely mediated through Af 3 - and A\j in Eqs.(1.2) and (1.3). 
In ID case, the complex coordinate Z{ is reduced to an ID variable, and it is easy to show 
that Afj and Afj vanish in Eqs.(1.4) and (1.5). Thus spinon and holon are decoupled in ID, 
and behave just like free-particles on their own tight-binding lattices. In 2D case, however, 
A\, and A\- can no longer be gauged away. A\- in Eq.(1.4) describes fictitious 7r-flux quanta 
attached to holons which are seen only by spinous in H s . Hence A^- will play the role to 
introduce doping effect into spin degree of freedom. On the other hand, A^- represents 7r-flux 
quanta bound to spinons which can be only seen by holons. A\- will then play the role of a 
scattering source in holon transport. So in the spin-charge separation scheme, new types of 
scatterings are present, and we will show in this paper that it is due to these unconventional 
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forces that magnetic and transport anomalies are produced in consistence with those found 
in the high-T c cuprates. 

Finally, let us outline how an electron is composed of the elementary excitations, holon 
and spinon, in this scheme. For both ID and 2D, the electron annihilation operator Cj CT will 
be rewritten as 



(-in 



hfho- 



;i.6) 



where 9i(l) is defined in Eqs.(2.19) and (2.20). Decomposition (1.6) can be understood as 
that annihilating an electron is equivalent to creating a holon and annihilating a spinon, 
and, at the same time, inducing overall nonlocal phase-shifts [in the brackets of Eq.(1.6)]. 
There are several ways to interpret the involvement of phase shift fields in Eq.(1.6). Since 
the nonlocal fields appearing in Eq.(1.6) resemble the Jordan- Wigner transformation in both 
ID and 2D cases,!!!! one may interpret them as statistical transmutations. Here hf is a 
hard-core bosonic operator, and bi a also satisfies bosonic commutation relation for the same 
spin index a but anticommutation relation for opposite spins [for details see Sec. II]. Thus 
these nonlocal fields in Eq.(1.6) are to guarantee Cj CT obeying an electronic commutation 
relation. A deeper physics behind this is related to the phase-shift idea, whose important 
role in a strongly-correlated system has been first realized by Anderson.il0@ In ID it has 
been explicitly shownll that an overall adjustment of the system occurs when a hole is 
doped into a large-?/ Hubbard chain, in order to retain the Marshall sign ruleS in the spin 
degree of freedom which is decoupled from the charge degree of freedom. This adjustment 
is found to be just represented by a phase-shift shown in Eq.(1.6), and it can lead to the 
correct Luttinger liquid behavior of the single electron Green's function. Furthermore, a 
phase-shift in ID can also be interpreted as reflecting the fact that each holon carries a 



spin domain wall [see Ref. [f2| and Sec. III]. Here one important distinction from Anderson's 
original phase-shift idea is that the many-body phase-shift in ID will not only give rise to the 
right non-Fermi liquid behavior, but also shift the Fermi surface to the position satisfying 
the Luttinger theorem. In other words, there is no need to assume a right Fermi surface at 
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the very beginning. In the present scheme, the form of Eq.(1.6) can ensure a large electron 
Fermi surface for both ID and 2D. But we shall discuss these electronic properties in a 
separated paper for the sake of clarity. In that same paper, we also show that when holons 
and spinons are both Bose-condensed, the phase-shift field in Eq.(1.6) will lead to a long- 
range pairing correlation of electrons (i.e., off-diagonal long-range order), and determine the 
symmetry of the superconducting order parameter. In the present paper, we mainly focus 
on the spin-charge separation formalism and explore the corresponding spin dynamics and 
transport properties. 

The remainder of the paper will be organized as the following. In section II, we construct 
the above-described spin-charge separation scheme!! based on the t — J model. At half- 
filling limit, a long-range AF order can be recovered in 2D, while the well-known spin-spin 
correlation is also obtained in ID. At finite doping, the gauge fluctuation is shown to be 
suppressed so that one has a real spin-charge separation. In section III, the spin dynamics in 
this scheme is studied at finite doping. In ID, the correct spin-correlation is to be reproduced. 
In 2D, a number of interesting properties are discovered which provide a consistent picture 
for the anomalous spin properties in the cuprates. In section IV, the transport properties 
is investigated. We show that there exists an interesting scattering mechanism within the 
present framework which leads an effective long-wavelength, low-energy Lagrangian. Such 
a Lagrangian can give rise to a series of exotic transport properties, in excellent agreement 
with the experimental measurements. Finally, a conclusive summery will be given in section 
V. 

II. SPIN-CHARGE SEPARATION: FORMALISM 

In this section, a mathematical formalism of the spin-charge separation state will be 
constructed based on the t — J model. 
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A. The slave-boson formalism and mean-field states 

According to Anderson,!! Zhang and Rice,@ low-energy physics in the cuprates layers 
may be properly described by a single band t — J model. The t-J Hamiltonian readsil 

H t _j = H t + Hj, (2.1) 

with a hopping term 

H t = -t £ c+c ja + H.c, (2.2) 

<ij><J 

and a superexchange term 

Hj = J £ (S* • S i - ^f) , (2.3) 

<ij> V 7 

in which S, = | X)<j ctrityva'Cia' ( 0" is the Pauli matrix), and rii = J2a c tcr c iv The Hilbert 
space of the Hamiltonian Eq.(2.1) is restricted by the no-double-occupancy constraint 

E^<1, (2-4) 

(7 

which imposes a strong correlation on electrons. 

By using the slave-boson decomposition of the electron operatorll'0 

c ia = fi a , (2-5) 

the no-double-occupancy constraint (2.4) can be replaced by an equality constraint 

hth l + J2f?Ji« = h (2.6) 

a 

where hf is a bosonic creation operator and fi a is a fermionic annihilation operator, known 
as holon and spinon operators, respectively. In this formulation, the hopping and exchange 
terms of the t — J Hamiltonian become 

H t = -t £ h+h i f£f jo + H.c, (2.7) 

<ij>cr 

and 
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Hj - E fiafitrfpfitT'- (2-8) 

<ij><J,cr' 



In obtaining Eq.(2.8), a term oc J5 2 (S is the doping concentration ) has been neglecte 
as usual for simplicity. 

The advantage of using the slave-particle formulation for the t — J model is that once the 
initial state satisfies the constraint (2.6), the no-double-occupancy condition can be always 
preserved under the Hamiltonian (2.7) and (2.8). Thus the constraint (2.6) in the slave- 
particle representation will not play a crucial role as the constraint (2.4) does in the original 
Hamiltonians (2.2) and (2.3). Furthermore, in the slave-boson formulation the natural mean- 
field decouplings usually have simpler structure at finite doping than those in the so-called 
slave-fermion, Schwinger-boson formalism.!! 

Various RVB-type mean-field states have been obtained in the slave-boson formulation.ill 
The mean-field decoupling of the Hamiltonians (2.7) and (2.8) may be realized by introducing 
the following mean-fields 

x£ =< f$,U° >; (2-9) 

H i:j =< hfhj > . (2.10) 

By neglecting higher order fluctuations around xfj an d Hij, the mean-field versions of Ht-j 
can be obtained in the form of H^ij = H^ IF + H^ IF , where 



H A h IF = -t £ x'ihthj + H.c- (2.11) 

<ij>a 



H? IF = ~^ E (xi + jH^ftJ^ + H.c. (2.12) 

<ij>acr' 

The constraint (2.6) should be simultaneously relaxed at this level to a mean-field one 

< hfhi >+Ea< ftJic >= L 

As shown by Eqs.(2.11) and (2.12), both H^ tF and H^ IF have similar hopping forms 
decided basically by one mean- field xlj (usually oc Xij a ^ the saddle-point). Since 
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Xij is responsible for both the hopping and antiferromagnetic superexchange strengths in 
Eqs.(2.11) and (2.12), an optimization of these two competing charge and spin processes 
may be relatively easy within this framework. 

A simplest mean-field state can be obtained by choosing xij = Xo and Hij = Hq, known 
as the uniform RVB state.il Such a state recently has attracted intensive attention, as one 
has been able to go beyond the mean-field (or saddle-point) approximation by including 
the phase fluctuations of xij an d in terms of the gauge-field descriptionliHi] A lesson 
learned from the gauge theory approach is that a singular transverse gauge fluctuation could 
become very important, which implies that spinon and ho Ion are usually confinedi! by gauge 
field, in opposing to the artificial decomposition of electron operator in scheme (2.5). The 
non-conventional transport phenomena have been discussedES in this framework which are 
both theoretically and experimentally interesting. 

The uniform RVB state is presumably energetically favorable at a sufficiently large dop- 
ing. On the other hand, when it is close to half-filling and the doping concentration S is low, 
the uniform RVB state may not be an appropriate saddle-point state because the antifer- 
romagnetic energy could be underestimated. For example, at half-filling limit with 5^0, 
the variational superexchange energy of the uniform RVB state, which can be computed by 
using the variational Monte Carlo ( VMC ) method with the no-double-occupancy condition 
(2.6) being exactly implemented, is found^l to be —0.53 J (per site), about 20% higher than 
the best ground state energy —0.6692 J [Ref. |B"4|| . Furthermore, no AF long-range order 
would appear in this limit as it should be present in the true ground state of H t _j. 

A substantial improvement of the variational ground-state energy at half-filling can be 
realized by introducing a phase to xij m Eq.(2.9), i.e., 

4 = Xoe* y . (2.13) 

(Note that = at 5 = 0.) Here the phase 9^ can not be simply reduced to a difference 
like 9i — 9j, or in other words, a summation of 9ij along the closed bonds of a plaquette, 
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is none zero. $□ may be regarded as some fictitious flux threading through the plaquette. 
At half-filling, the variational energy is optimized at $ n = 7r, known as the 7r-flux phase,H 
whose VMC value of the exchange energy is —0.623 J [Ref. [63|], only 5% higher than the 
exact ground-state energy. However, the absence of a long-range AF order in the 7r-flux 
phase implies that the long-wavelength, low-energy AF correlations are still underestimated 
which should predominantly account for the 5% energy missing at this saddle-point state. On 
the other hand, the high-energy properties of the Heisenberg model, like the temperature 
dependence of the specific heat may be well explained in terms of the 7r-flux phase, as 
discussed in Ref. |65| (In Ref. |65], the 7r-flux phase is obtained from a different formulation, 
involving a 2D Jordan- Wigner transformation^! of the spin-1 operators, instead of the spin-^ 
operators used in the present framework.) 

Thus, by comparing with the uniform RVB state at half-filling, the short-range AF cor- 
relation in the 7r-ffux phase has been substantially improved as a result that the transverse 
gauge fluctuation in the former is condensed as the static uniform flux in the latter case. 
Nevertheless, lack of long-range AF order in the latter means that the long-range AF cor- 
relation is still not appropriately accounted here. In the following, we will consider a new 
type of saddle-point states which could retain the high-energy characteristics of the it— flux 
phase, while properly include the long-wavelength AF correlation in low-energy regime. 



B. Flux binding: new saddle-point 

In the 7r-flux phase discussed above, the mean-field Hamiltonian describes a non- 
interacting fermion gas under some uniform fictitious magnetic field. However, according 
to the exact diagonalizations on small latticeJH such a type of states is generally not en- 
ergetically favorable as compared to a state in which the uniformly-distributed lattice flux 
is quantized into infinitesimal-size flux tubes, and each of them is bound to an individual 
particle. The former may be regarded as a "mean-field" version of the latter. In general, 
the latter is known as an anyon system^ as the new composites of particles and flux-tubes 
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obey different statistics, depending on the ffux strength of each flux-tube, as compared to 
the statistics of the underline particles. 

This numerical result is very suggestive here, as in the slave-boson formalism there ex- 
ists a gauge (phase) degree of freedom which can allow the above flux binding procedure 
to happen if it is indeed energetically favorable. This gauge freedom is manifested in the 
decomposition (2.5), where one may always associate phases e ld and e~ %e to hf and fi a , 
respectively, without changing Ci a . Such a freedom reflects the arbitrariness of the decom- 
position (2.5). A gauge field will then play a role to confine any non-physical spinon and 
holon together as an electron. Only when one finds a correct spin-charge decomposition, 
the gauge fluctuation can be expected to get suppressed. Our strategy in the following is 
to optimize the flux phase by exploiting such a gauge freedom. Main procedure is to regard 
the uniform fictitious flux as a "mean-field" version of the system where fluxes are quantized 
and bound to the particles. This latter system may be generally called as a flux-binding 
state. Given the above-discussed uniform 7r-flux, there can be two ways to construct the 
corresponding flux-binding states at half-filling, as will be outlined below. 

Scheme one All of no matter their spins, are bound to the same type of flux 
quanta. When these flux quanta are uniformly smeared out in space, one should recover the 
uniform flux in 7r-flux phase. At half-filling, the total number of spinons is one per plaquette 
on average. Thus each flux tube has to be quantized at 7T, such that 

lea <j 

with n{ a = ffofia- On average < J2n @ij >— 11 ■ 111 this scheme, since each flux-tube is 
quantized at tt, spinons are effectively turned into semions0 (an exchange of two semions 
gives rise to a phase Scheme one has been already discussed in detail in Ref. Away 
from half-filling, the mean-field version of this flux-binding state corresponds to the so-called 
commensurate flux-phase (CFP) [Ref. |6"5fl . And the gauge theoriesE^Hll based on CFP has a 
close connection with such a flux-binding state. In both cases, superconducting condensation 
is found in the ground state, and is related to the semionic condensation.0 It has been shown 
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in the flux-binding stated that there is no explicit time-reversal and parity symmetry broken 
due to the cancellation between charge and spin degrees of freedom. The most interesting 
features show up in the normal state,0 which can well explain the anomalous in-plane 
transport properties in cuprates, including resistivity, the Hall effect, and thermopower. 
But such a state still has trouble to recover the AF ordering at half-filling. 

Scheme two In this scheme, the spinons with opposite spins may not see the flux- 
tubes carried by each other. Thus the flux quantum of each flux-tube has to be doubled in 
order to recover the 7r-flux phase in the mean- field limit. In this case, the phase 9ij as seen 
by the spinons with spin index a has to satisfy the following condition 

E D ^ = 27r E^L (2-16) 

lea 

such that < J2n^ij >= 71 & t half-filling. Equation (2.16) means that each spinon is bound 
to a 2tt flux-tube, which would effectively change the spinon statistics from fermion into 
(hard-core) boson, for spinons with the same spin index a. In this second flux-binding 
scheme,il the experimental features of the transport properties found in Scheme one will be 
essentially retained, while interesting magnetic properties can be also obtained here. As an 
example, the AF long-range order can be recovered at half-filling. Furthermore, the correct 
ID behaviors is to be reproduced naturally within this framework. All of these suggest 
that Scheme two may be a better approximation than Scheme one for the same physical 
state. Therefore, we shall fully focus on the saddle-point of Scheme two later, and explore 
its various properties in the remainder of the paper. 

1. Mathematical scheme: Half-filling 
At half-filling, H t has no contribution. One may rewrite Hj in Eq.(2.8) in a form 

Hj = ~ E (e^f^f ja )(e<f^f lP ), (2.17) 

Z <ij>a/3 V J 

where a phase A?j, which will play a role similar to in Eq.(2.16), is introduced: 
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which depends on the positions of spinons nonlocally. This nonlocality is introduced through 
the multiple- valued phase 

9 i (l) = Im\n(z i -z l ), (2.19) 

with Zi = Xi + iyi in 2D. A vorticity ±2n can be obtained in if the coordinate Z{ has 
been continuously changed around z\ once. In ID case, 6i{l) in Eq.(2.19) will reduce to 

6i(l) = ±n6(l-i), (2.20) 

where 6(x) on the right- hand- side is the step function. In obtaining Eq.(2.17), the constraint 
condition (2.6) has been used, i.e., 

E n £ = 1 » ( 2 - 21 ) 

a 

at 5 = 0. 

The identical transformation in Eq.(2.17) has no real physical meaning until the mean- 
field decoupling scheme is introduced. One may define a hopping operator 

x% <' v >f,-fj»- (2.22) 

Then under the mean- field Xo =< Xij >> a saddle-point Hamiltonian of Eq.(2.17) is found 

as 

H s = -J X o E eW&fi* + H - c - (2-23) 

<ij><7 

Now the nonlocal field appears as a phase in the hopping matrix, which is equivalent to 
the role of ^ in Sec. A 

Two-dimensional case In terms of Eqs.(2.18) and (2.19), a summation of the phase 
A?j along a plaquette gives 

E D A* = 2naY.nl - 2n5 a , v (2.24) 
lea 
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Besides a (— 27r) lattice flux for a =|, which has no real physical meaning, the right-hand- 
side of Eq.(2.24) implies that each spinon of spin a carries a flux-tube quantized at 2tt, with 
a sign a — ±1. Since < n la >= 1/2, one has J2n < >— — 71 ■ That is, the "mean- 
field " version of the saddle-point Hamiltonian (2.23) just corresponds to a 7r-flux phase 
Hamiltonian as required (note that — ir and it fluxes per plaquette are equivalent). 

As pointed out before, the flux binding procedure discussed here is related to a statistics- 
transmutation which is easy to see after introducing a new operator bi a 

K = f la e- l ^ m ^-^ ] . (2.25) 

Then H s is reduced to 

H s = -J X o E (-<r)b&h* + H - c -i (2-26) 

<ij>a 

One can easily check that bi a satisfies the hard-core boson commutation relations (for the 
same index a ): \b ia , bf a ] — (i j), etc. Of course, for opposite spins one still finds 
anti-commutation relations. 

Equation (2.25) resembles the Jordan- Wigner transformation in both ID and 2D 
casesSil which changes the statistics of a fermion into that of a hard-core boson. Thus 
the presence of a uniform lattice flux in the 7r-flux phase is a precursor for a fermionic 
spinon to become a boson. Such a "bosonization" tendency for spinons in flux phase has 
been already noted before.^! Bosonic representation has been known generally to be superior 
in the description of spin AF correlations. 

Equation (2.26) can be further written in the form of a standard hopping Hamiltonian 

H s = -Jxo E bfb ja + H.c. (2.27) 

<ij>u 

by redefining 

b la = (-c?) % b ia , (2.28) 

in which o = ±1, and (— l) 1 = +1 for the even-sublattice and —1 for the odd-sublattice of 
a bipartite lattice. The spin operator like = f^fii can then expressed in terms bi a as 
follows 
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S+ = (-l)%bn. (2.29) 

One also has = (St ) + , and Sf = \Y,o °~Ko-bio- 

The Bose-condensation of bi a at zero-temperature as determined by Eq.(2.27) means 

< St >= (-IT < b+ >< b n 0. (2.30) 

In other words, an antiferromagnetic long-range order in the x-y plane can be indeed obtained 
in the present flux-binding scheme. The direction of the magnetization in x-y plane will be 
determined by the relative phase between the up and down spinons in Eq.(2.30), and could be 
arbitrary. The evaluation of the magnetization value in Eq.(2.30) needs a detailed knowledge 
of the hard-core boson behavior governed by H s . Since we do not know the exact solution 
of the hard-core boson Hamiltonian (2.27), we can not directly estimate the variational 
ground state energy of the Heisenberg Hamiltonian Hj in the present saddle-point state. 
Nevertheless, the ground state energy for a hard-core boson system like Eq.(2.27) is generally 
lower as compared to the uniform 7r-flux system in the fermion representation, which has 
been discussed by Gros et al. in Ref. |7B[ As pointed out by these authors, the ground state 
energy of a hard-core boson system described by H s could be as lower as 11% than that in 
the 7r -flux-phase. Of course, one still needs to check the variational energy of the original 
Hamiltonian Hj in the present saddle-point state. By using a trial bosonic wavefunction for 



H s given in Ref. [76], which is obtained by taking the absolute value of the known ferminic 
7r-flux-phase wave function and gives 7.4% higher variational energy than the best estimate 
for H s , we found analytically the variational energy of Hj is just identical to that in the 
7r-flux phase. The latter is already known from the VMC calculation, i.e., 5% higher than 
the exact value of the Heisenberg Hamiltonian. Since the above-mentioned bosonic trial 
wavefunction does not have a long-range order (Bose-condensation), one may attribute the 
5% higher ground state energy to it. A better bosonic trial wavefunction of H s should further 
improve the variational ground state energy of Hj. 

We note that beyond the present mean-field approximation, a gauge fluctuation should 
be considered. As will be discussed in Sec. C, such a gauge fluctuation is not gapped at 
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half-filling and thus could become very important. In fact, spin- 1/2 spinons will presumably 
be confined by the gauge field to form spin-1 excitation (spin wave) at half-filling. This 
bosonic description is different from the Schwinger boson representation. The latter has 
been found^ to well describe the low-lying spin excitation in 2D half-filling, even at a mean- 
field level. In contrast, the present scheme will become more powerful in spinon deconfined 
cases. For instance, it will be shown later that the gauge fluctuation is to be suppressed 
at finite doping, and with the presence of spin-charge separation there, the present bosonic 
representation can become quite convenient to include doping effect. An another spinon 
deconfinemnet case is in ID, which will be discussed below. 

One dimensional case As a further check of the present saddle-point state, we now 
turn to the one dimensional case where the asymptotic spin-spin correlation of the Heisenberg 
chain has been known for many years.@ 

Under the saddle-point Hamiltonian (2.27), one can write 

< S+(t)S-(0) >= (-1)'-' < b+(t)b n (0) >< &u(f)&£(0) > • (2-31) 

In order to evaluate the averages on the right-hand-side, the behavior of the hard-core boson 
described by foj(t) has to be determined first. A trick which has been often used in ID is 
to transform the bosonic operator b{ into fermionic operator again. By using the expression 
(2.20) for 6i(l) in ID, 6 i(T is expressed in terms of fi a in Eqs.(2.25) and (2.28) as follows 

b ia = fi^ ia ^>^. (2.32) 

One also finds that Afj = in Eq.(2.23), and thus in the fermionic representation H s simply 
describes a free lattice fermion gas. Then the correlation function like < b%(t)bj^(0) > in 
Eq.(2.31) can be straightforwardly calculated in long-time and -distance as follows^ 

< 6}(t)6 it (0) >= < /+(t) e ± -S i > I «f T W e ^E i > J »f T (o) /jT(0) > 

1 

°" [(Xi-Xtf-Vji 2 ]!/*' (2 ' 33) 

with Vf = Jxo a - The average < 6jj_(t)6^_(0) > shows the same asymptotic form. Finally by 
using (— l)' t_J = cos^(xi — Xj) and the rotational invariance, one obtains 
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which describes the correct asymptotic spin-spin correlation in ID and was first derived by 
Luther and Peschel from XXZ model.0 

Thus the present flux-binding saddle-point state has produced a correct result in ID, 
even though it is originally constructed in 2D. This is in contrast to the Schwinger-boson 
approach,!! which works rather successfully for the 2D Heisenberg model, but could not 
predict the correct behavior for the ID spin 1/2 case. This may be attributed to the fact 
that spin- 1/2 excitations exist in ID Heisenberg chain, but not in 2D. The present bosonic 
representation properly describes these spin- 1/2 spinons. 

We note that in the present approach the spin-spin correlation in z-axis can not be 
evaluated at the same level of approximation as in the x — y plane. The reason is that 
< Sf ■ Sj >~ | J2a < nb icr nb j<T > basically involves the density- density correlation of spinons 
with the same spin index, while < Sf ■ Sj > is related to the single-particle propagators 
of spinons under the saddle-point (2.27). The former is generally much harder to evaluate 
because of the hard-core condition of bi a . The ID correlation < Sf ■ 5| >~ (— /\xi —Xj\ 
could be obtained with the hard-core condition being accurately treated through the Jordan- 
Wigner transformation. 

One final point we like to stress here is that H s in Eq.(2.23) is formally the same as that 
in the uniform RVB saddle-point, since A\- = in ID. In other words, the flux-binding, or 
the statistics-transmutation, has no effect on the Hamiltonian in ID as is well known. What 
is different between these two saddle-point states is the boundary condition for the fermion 
f ia : in the present saddle-point, f ia appearing in Eq.(2.32) will become multiple- valued at 
the boundary once the ID chain is closed as a loop. As discussed by Shastryjll such a 
boundary condition decides why a hard-core boson system is energetically more favorable 
than a fermionic gas. 
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2. Flux-binding: finite doping 



For a small doping concentration, one expects the antiferromagnetic correlation still to 
remain dominant, even though it may be strongly modified by doped holes. It is very 
interesting to see the evolution of the flux-binding saddle-point state in this finite doping 
case. 

We may follow the procedure shown in Sec. 1, starting with the identical transformation 
of Hj in Eq.(2.17). At finite doping, the form of Afj in Eq.(2.18) will be slightly modified 
because the constraint Eq.(2.21) now is replaced by the full equation (2.6). It is easy to 
show that A\a is given as follows: 




(2.35) 



which keeps the identical transformation in Eq.(2.17) unchanged. In Eq.(2.35), an additional 
factor nf/2 is introduced, with n\ = hfhi as the ho Ion number operator. 

At finite doping, the hopping term H t in the H t _j will also contribute, and it may be 
rewritten in a form 

H t = -t £ e<K} h,\h- (2-36) 

<ij>cr 

where yf-- is defined by Eq.(2.22) with A°- given in Eq.(2.35). A{j in Eq.(2.36) is a re-arranged 
form of A\j after using the constraint (2.6): 

4 = \H P*(0 - d M fe ™l ~ • ( 2 - 37 ) 

Note that A{j in Eq.(2.37) has no explicit a-dependence. 

Now we may take the mean-field decoupling of Hj in Eq.(2.17) and H t in Eq.(2.36), by 
introducing the mean-fields Xo =< Xij > an( i Hq =< Hij > (here Hij = e tA ^hfhj). At the 
same time, the constraint (2.6) has to be relaxed up to a global level, as outlined before. 
Then a saddle-point Hamiltonian at finite doping is found to be 

= H S + H h , (2.38) 
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in which the spinon degree of freedom has the same form as in the half-filling case: 

H s = -J s J2 e^if+fjr + H.c, (2.39) 

<ij>,cr 

with J s = Jxo + tH . The holon degree of freedom is described by: 

H h = -t h e^h+hj + H.c, (2.40) 

<ij> 

with th = 2txo- th and J s have to be determined self-consistently. Since the averages xo an d 
H actually do not depend on t h and J s , the self-consistency should be always satisfied here, 
even though an actual determination of the values for th and J s is nontrivial. Generally one 
can estimate th ~ t and J s ~ J at small doping. In the bosonic representation of Eqs.(2.25) 
and (2.28), H s can be further expressed as 

H s = -J s e itTA *ibib ja + H.c, (2.41) 

<ij>(T 

where the topological phase is defined by 

4 = ^£^(0-^(0K- (2-42) 

In one dimension, one will find = A{j = in Eqs. (2.40) and (2.41) according to 
Eq.(2.20). The holon (described by hi) and spin (described by b ia ) degrees of freedom are 
thus decoupled at this saddle-point state. In two dimensions, however, the holon and spinon 
are coupled with each other through the nonlocal fields A{j and A^ in H h and H s . The 
topological gauge field A\- can be interpreted as describing fictitious 7r flux quanta attached 
to holons, but only seen by spinons. This becomes clear if one considers an arbitrary loop 
C, and counts the flux enclosed, which is given by 

Ec^^E^- (2-43) 
A{j in Eq.(2.37) may be rewritten as A{j = Afj — 4>%, where 

^ = lT,m)-dj(l)](T,^iX (2-44) 
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and 



is denned by 



4 = m)-om 



(2.45) 



which describes a lattice 7r-flux with J2n 4>% = n - The topological gauge filed A\- describes 
fictitious 7i flux-tubes carried by spinons which are seen by holons. In Sees. Ill and IV, we 
will investigate how the couplings induced by nonlocal fields A\j and Af- will lead to highly 
nontrivial spin and charge properties in 2D case. 

To end this section, we would like to give a different perspective about the construction 
of the flux-binding state. Recall that one has a gauge degree of freedom in a slave-boson 
decomposition: Ci a = h~tb ia = (hte l6h ^j f/io-e^-M, where Oh and Of can be anaphase satisfying 
Oh + Of = 0. Particularly, if Oh and Of are chosen as the topological phases in the Jordan- 
Wigner transformation, statistics-transmutation can happen so that one may end up with a 
slave-fermion, or more generally, slave-anyon decomposition.^ This is not surprising because 
all of these formalisms are mathematically equivalent. Much complicate decomposition can 
be also constructed. Distinctive physics is involved here only when one makes the mean- 
field decoupling. In this way, different decomposition means different saddle-point state. The 
usual slave-boson and slave-fermion formalisms are the simplest ones, but not necessarily the 
appropriate ones for the correct saddle-point. The present saddle-point state corresponds 
to Oh = l/2J2i^i0i(l)[(r(l — nf) + J2a an i a ] an d Of = —io' 0i(l)n{ a . It is easy to check 
that Oh + Of = in terms of the constraint (2.6). Correspondingly, Cj CT may be rewritten 
as in Eq.(1.6) where bi a is defined by Eqs.(2.25) and (2.28). In this new decomposition, an 
electron is composed of two bosonic holon and spinon, together with some nonlocal fields. 
The spin operator like = fXfi\ can also be expressed as 



Its physical meaning will be explored in Sec. III. In the following, we will go beyond the 
mean-field approximation and show that a real spin-charge separation can be indeed realized 
with the decomposition (1.6). 



(2.46) 
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C. Gauge-theory description: Spin-charge separation 

It is now a well-recognized factES that any saddle-point-state properties of the t — J 
model could be substantially modified by the low-lying gauge fluctuation around the saddle- 
point. So it is very important to check the gauge fluctuations around the present saddle- 
point. Here gauge fluctuations are to be related to phase fluctuations of the mean- fields 
< Xij > and < Hij >. The main idea behind the construction of flux-binding saddle-point 
states is to incorporate the most singular transverse gauge fluctuations into the saddle-point 
via flux binding so that additional phase fluctuations become insignificant. These singular 
gauge fields are represented by the topological gauge fields and in the saddle-point 
Hamiltonians H s and Hh- In the following, we will show that the transverse gauge fluctuation 
around this saddle-point is indeed suppressed at finite doping, and therefore we have a real 
spin-charge separation whose low-energy physics will be determined in terms of the effective 
Hamiltonians H s and Hh, defined by Eqs.(l.l)-(1.5) in the Introduction. 

An effective Lagrangian, with the constraint (2.6) enforced and the (phase) gauge fluc- 
tuation around the present flux-binding saddle-point state included, can be written down 
aSH 

£*ff = E (E b t(dr - Ah* + hf (d T - fj,)h^ + H eff , (2.47) 

where 

H*ff = -J. E e^+^bib^ + cc. 

<ij>ir 

-t h e^^-^hfhj + c.c. (2.48) 

<ij> 

Here in Eq.(2.47) describe an internal gauge field. Originally the no-double-occupancy 
constraint Eq.(2.6) [with changed to bi a by transformations Eqs. (2.25) and (2.28)] is 
implemented through a Lagrangian multiplier field, whose fluctuating part is then absorbed 
by the longitudinal part of in Eq.(2.48) due to the gauge invariance,il with \i left in 
Eq.(2.47) enforcing the constraint at a global level. 
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The existence of a gauge freedom in C e f / will guarantee the following current constraint 
between spinon and holon:iM! 

J s = -J h - (2.49) 

It is noted that Eq.(2.49) in longitudinal channel is a simple reflection of the density con- 
straint (2.6), but its transverse channel has nothing directly to do with the no-double- 
occupancy constraint, and is related to the property of the t — J Hamiltonian. It is the 
transverse gauge fluctuation that becomes singular in long-wavelength, low-energy regime 
in the uniform RVB state. This transverse field would serve as a confine force in a usual 
non-spin-charge-separation state. The current constraint (2.49) is also connected to the 
Ioffe-Larkin combination rulel! of the response to an external electromagnetic field 

K e = [K; 1 + K^]~\ (2.50) 

where K s and are the response matrices of spinon and holon systems under effective 
Hamiltonians H s and in Eqs.(1.2) and (1.3). 

In gauge theory, the dynamics of the gauge field is determined after the spinon and 
holon degrees of freedom are integrated out .00 When the gauge fluctuation is weak, one 
may use the Gaussian approximation, and the gauge field propagator D a = — < T T a^a v > 
(//, v — x, y) in imaginary-time can be determined by@@ 

D a = -[K s + K h ]-\ (2.51) 

in imaginary-frequency space. In the present system, holons are under some fluctuating flux 
described by Afj [cf. Sec. IV], and its response function will follow a usual metallic q and u 
dependances:0 ico&h — XhQ 2 - The most interesting behavior will come from the spinon part 
K s as discussed below. 

In the spinon part of H e ff, there is a sign a = ±1 in front of the topological phase A\y 
This sign means that spinons with \ and \ spins see the fictitious flux quanta carried by 
holons in opposite directions. Since there is no such a sign in front of the gauge field a^ , a 

27 



non-zero will then polarize the spinon system with regard to spin a. This polarization 
is generally energetically unfavorable, and we expect a suppression of the gauge field to 
stabilize the system. 

A mathematical demonstration is straightforward. As will be discussed in Sec. Ill, 
one may rewrite A^- as A-,- = A 1 }- + 5A^-, where A- - is a "mean- field" one with the flux- 
quanta uniformly smeared out in space. A uniform fictitious magnetic field Bh = vr5/ a 2 will 
correspond to such a vector potential A^. On the other hand, bA\- will be correlated with 
the holon density fluctuation. If one integrates out the spinon degree of freedom (under the 
"mean-field" A^) in £ e //, an effective action for quadratic fluctuation of the gauge field is 
found by 

S s [a, 5A h ] = - U + ^A h ) n ff (a + o5A h ) , (2.52) 

where il^ is the response matrix of spinon a under a fictitious magnetic field aB h . K s will 
be determined through S s = ^aK s a [Ref. [7C[| , and one needs to further integrate out 5A h in 
Eq.(2.52). At first step, if one neglects 5A hj s own dynamics, by directly integrating out SA h 
in Eq.(2.52), the following expression for K s can be obtained 

(2.53) 

The key thing here is that the up and down spinons see the opposite fictitious magnetic 
fields, and the Hall conductance terms in Ilf and Il| will have opposite signs. As shown by 
Wiegmann0 in a generic situation, due to the same amplitude but opposite signs of the Hall 
conductances, K s in Eq.(2.53) must exhibit the "Meissner effect" in the transverse channel: 
i.e., Kg = l/47rA 2 7^ at uo = 0, and q — > limit. Here we find 1/A 2 = nBhJ s - Such a 
"Meissner effect" implies that the gauge field will get suppressed in the spinon system. In the 
long- wavelength and low-energy limit, K s becomes dominant in the transverse channel of D a 
as Kh ~ 0, and a non-zero Kj- then leads to a gap oc 5 J in the transverse gauge fluctuation 
of Eq.(2.51) (Anderson-Higgs mechanism). This gap means a spin-charge deconfinement as 
pointed out before. In contrast, a finite K~ l will play a negligible role in the electromagnetic 
response function K e in Eq.(2.50) at small q and a>, and one finds 
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K e ~ K h , 



[2.54) 



in such a limit. So the holon degree of freedom solely determines the long- wavelength, low- 
energy response to an external electromagnetic field in this spin-charge separation system. 
Obviously, all these features will disappear at 5 — > as vanishes. 

The above conclusion is still true if one includes SA h, s own dynamics by adding a term 
\8A h D A 5A h in Eq.(2.52). Here D A is the free propagator for 5A h which is related to 



holon density fluctuation as follows: 
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-D h 

,2 C 1 



(2.55) 
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with D h as the holon density- density correlation function. Then Eq.(2.53) is modified as 
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Uf +Ul 1 (D A )- 1 U^ 1 



(2.56) 



Since D A does not have a transverse- longitudinal mixed term (the Hall term), the above 
conclusion about K s can be easily found to remain unchanged as long as D A approaches to a 
constant or vanishes slower than q 2 when q — > 0. (Note that the density-density correlation 
function^ of a hard-core boson gas is very similar to a free-fermion gas, and in the latter 



case one has D h ~ const, at u = and q — > 0.) 



III. SPIN DYNAMICS AT FINITE DOPING 

Spin dynamics in the insulating cuprates has been well understood within the 
framework0 of the Heisenberg model. A real important issue is how the antiferromag- 
net is affected under doping. In this section, we will explore such a doping effect under 
the spin-charge separation scheme, and compare its unique features with those found in the 
cuprates. We first consider the ID case, where the exact analytical results are available for 
comparison. 
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A. Spin-spin correlation in one dimension 



The present case corresponds to the large-?/ limit of the Hubbard model, whose ex- 
act solutionS has been obtained a long time ago. However, only recently the underlying 
physics and various correlation functions have been clarified, following Anderso rffland 
other authors.iB@i 

The spin-spin correlation function is one of important correlation functions in ID. Ac- 
cording to the numerical calculation^ based on the exact solution, a 2k j oscillation is present 
in the spin-spin correlation function. One may naively relate this incommensurate structure 
with the existence of a large electron Fermi surface (points at kf and —kf). Similar argument 
may also be applied to a 2D case. This observation should be physically reasonable, but by 
itself is not sufficient to get the correct spin dynamics. Much rich effects will be involved 
here due to strong electron-electron correlations. In fact, one has a spin-charge separation. 
In this case, the spinons are presumably responsible for the spin dynamics. But at the same 
time, holons as solitons of doped holes bound with spin domain-walls will also contribute 
to spin properties. A combination of these decides a peculiar doped spin-spin correlation, 
which is highly non-trivial from a Fermi-liquid point of view. In the following, we will show 
how a correct spin-spin correlation can be obtained without directly involving the electron 
Fermi-surface. 

In the present scheme, the spin flip operator Sf in Eq.(2.46) can be expressed in ID in 
terms of Eq.(2.20) as 

st = (-iyb+b lie ^>^ . (3.i) 

Here a prominent feature is a Jordan-Wigner-type nonlocal phase 7r J2i>i itf which involves 
holon number operator. At first glimpse, one may feel strange for a holon number operator 



to appear in a spin expression. But this is not new, and has already been found in Ref. [12 
by a path-integral approach, as a result of spin-charge separation. It actually describes the 
effect of spin domain-walls carried by holons mentioned above. This is easy to see if one 
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freezes the dynamics of spinon b ia by letting it "condensed" in Eq.(3.1). Then one would 
have a Neel order in spin x — y plane when holons are absent. For each added holon at I site, 
one finds an extra sign e ±l7T = — 1 for every at i < I, which means a flip for all those spins. 
In other words, there is indeed a spin domain-wall accompanying with each doped hole, and 
that is the physical reason why the holon number operator enters into Eq.(3.1) nonlocally. 
We note that the nonlocal field appearing in Eq.(3.1) can be also related to the phase-shift 
field that leads to the Luttinger liquid behavior and an electron Fermi-surface satisfying 
the Luttinger liquid theorem in the single electron propagator, which is to be discussed 
elsewhere. In fact, a 2kf oscillation in the spin correlation function will arise naturally due 
to this nonlocal field as shown below. 

Due to Afj = = in ID, holon and spinon are decoupled in the Hamiltonians (1.2) 
and (1.3). The nonlocal phase in Eq.(3.1) will then solely determine the doping effect. One 
may write 

< #(f)ST(0) >= (-1)*-' < &+ T WMo) >< &,-i(f)&5i(o) > 

x < e ± "E I >."fWe Ti ^ I >/"(°' > . (3.2) 

Since bi a and hi as described by Eqs.(1.2) and (1.3) are hard-core bosons, each average on 
the right-hand-side of Eq.(3.2) can be evaluated straightforwardly in the asymptotic limit 
(see Sec. II B, and Ref. j42|) . For example 



< e^E^nfWg^Ew^) >oc (3.3) 
with x = Xi — Xj, and Vh = 2i/ l asin(7r5). The final result is 

< Sf{t)ST{0) >0C (a;2 _ t; 2 t2) l/ ( 2(a f 2 l t; 2 t2) l/4' ( 3 - 4 ) 

where kf — (1 — 5)7r/2a. The doping effect as described by the phase-shift in Eq.(3.1) enters 
Eq.(3.4) by changing the oscillation from 27r/a to the incommensurate 2kf, and at the same 
time contributing an additional power (x 2 — t^t 2 )" 1 / 4 . It is noted that < Sf (t)S , |(0) > can 
not be directly computed at the same level of approximation as explained in Sec. II, but it 
should follows the same behavior due to the rotational invariance. 
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Equation (3.4) recovers the correct spin behavior at the strong-coupling fixed point of the 
model (i.e., the large-?/ Hubbard model). It suggests that the present state indeed catches 
the correct characteristics of charge-spin separation in ID, even though such a scheme is 
originally constructed for 2D. 

B. Doping effect in two dimensions 

In the one-dimensional case, spinon and holon are decoupled, and the doping effect on 
spin dynamics is solely contributed by the spin domain walls carried by holons. In 2D, this 
nonlocal doping effect takes a different form, because spinon and holon will no longer behave 
like free solitons. As noted before, the topological phase in H s cannot be gauged away 
in 2D, and it will represent a new type of nonlocal influence of doping on spinons. 

1. New length and energy scales introduced by doping 

As shown in Sec. II, the gauge field confining spinon and holon are suppressed at finite 
doping. But there are still residual interactions between spinons and holons, and spinons 
can always feel the existence of holons nonlocally by seeing the flux-quanta bound to the 
latter. 

Let us first see how this exotic interaction can change the topology of a holon. In the 
spin flip operator (Eq.(2.46)), one has a nonlocal phase [X^Im In (zi — z^n*)] involving 
holon number operator. In resemblance to the domain-wall picture in ID, one may interpret 
it as describing a spin vortex in x — y plane with vorticity 2n bound to each holon. This 
could be seen if one freezes the dynamics of bf a and treats it as a number in Eq.(2.46). 
But bf a here can no longer be regarded as a single-valued quantity because it is under the 
topological phase in H s . A Berry-phase counting shows that the phase in bf a cancels out 
the effect of the nonlocal phase-shift in Eq.(2.46) such that there is actually no 2n- vortex 
topological texture formed around a holon. Nevertheless, if one goes along a line across such 
a holon, one can still find a domain-wall-like spin singularity at holon site similar to that 
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in ID. Thus, a holon in 2D is no longer associated with a spin topological object. Instead, 
it may carry a Shraiman-Siggia-typelll dipolar texture with vorticity = 0. Such an object 
will not break the T- and P-symmetries. One expects a strong dynamical renormalization 
to be involved in determining the profile of each holon. However, we are not interested in 
a single-doped-hole problem here. We shall focus on a finite doping concentration in the 
following, where a useful mathematical description becomes available. 

For simplicity, we are going to use a continuum version of the spinon Hamiltonian H s by 
taking the lattice constant a — > 0. This continuum approximation may be justified at small 
doping and in the low-energy, long-wavelength regime, where the amplitude of is small 
and the spinons as bosons mainly stay near the bottom of energy band. Such a continuum 
version for H s in Eq.(1.2) can be easily written down as [cf. Ref. Q 

H s = E/ d 2 r K(*) HV ~ m a s Ah)2 K(r), (3.5) 

with m s = (2J s a 2 ) -1 and A h (r) defined by 

1 r ,o ,z x fr — r' 



A ft r 



1 r , 9 ,z x (r — r , ,„ . 

2j dr | r _ r 'P Ph{r) - (3 ' 6) 



Here p^(r) is the holon density ph = h + (r)h(r). 

We note that Eqs.(3.5) and (3.6) would describe a semion problem^ if the holon density 
Pft(r') were replaced by the spinon density. The mathematical similarity suggests that one 
may borrow the method developed in the anyon problem.fiiilH The idea is to rewrite = 
Ph + (Ph — Ph) m Eq.(3.6) by introducing an average holon density p\ = 5 /a 2 . Then A h (r) 
is rewritten as 

A h (r) = A h (r)+5A h (r) . (3.7) 

Here A h (r) corresponds to p^, and in the symmetric gauge it may be expressed bySii 
A h (r) = f^(z x r) with Bh = irph- Physically, the vector potential A h describes a mean- 
field magnetic field Bh obtained after the flux quanta bound to holons are uniformly smeared 
out in space. SA h (r) in Eq.(3.7) is the fluctuation part 
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6A h (r) = \ J d ^^ r _^\ Ph {v>) - p h ), (3.8) 

which can be treated perturbatively. We emphasize that in an anyon problem, 5A h would 
correlate with anyon density fluctuation and representSii a long-range interaction among 
anyons. Here 5A h is determined by the density fluctuation of holons instead of spinons 
themselves, which belong to an independent degree of freedom. Thus one expects 5A h (r) to 
provide an independent dynamical scattering source in Eq.(3.5) just like phonon in a usual 
electron system. 

The separation (3.7) is meaningful when the hole density is not too low. With the 
presence of a fictitious magnetic field Bh, a new length scale is introduced to the spinon 
system, which is the magnetic cyclotron length 

l c will later be connected to the antiferromagnetic correlation length. At this mean-field 
level, a Landau-level structure appears in the spinon energy spectrum, and a basic energy 
scale is the cyclotron energy 

u c = — = 2n5J s . (3.10) 

Another energy scale measuring the broadening T s of each Landau level caused by the 
fluctuation 5A h should be also correlated with the doping concentration. T s will be related 
to an important low-energy scale in spin dynamics. It can be estimated to be in an order 
of magnitude ~ 5 J if the energy scale of the fluctuating 5A h is sufficiently small. T s could 
be even sharper when the fluctuating 5A h has a higher energy scale as holons become more 
mobile at finite doping. Of course, different Landau levels generally could have different 
broadening widths, and a further discussion will be given later. Similar broadening problem 
in a semion system has been recently studied by Levy and Laughlin,0 where the dynamics of 
5A h is already known from the RPA calculation. In the present case, the dynamics of 5A h is 
directly related to that of 5ph = Ph~ Phi which will be in turn determined by coupling with 
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the spinon degree of freedom. So a self-consistent treatment of the Landau level broadening 
could be much complicated here. Nevertheless, for the purpose of understanding basic 
characteristics for spin dynamics, only a general shape of the Landau level broadening is 
needed, and the detailed structure will not be crucial. 

Besides the length and energy scales, we point out that there also exists a basic temper- 
ature scale which is related to the Bose condensation temperature T* of the bosonic spinon 
b ia . Recall that for a 2D free boson gas, a 2D Bose condensation is always suppressed by the 
thermal excitations at any finite temperature, due to finite density of states at low energy. 
In the present case, due to the broadening of the lowest-Landau-level (LLL), the low-energy 
density of states presumably will fall off continuously to zero at the LLL energy bottom. 
Thus the low-lying thermal excitation has vanishing weight at the low-energy tail, which 
could not kill the Bose-condensation of spinons at a sufficient low temperature. As will be 
demonstrated later, T* will represent an important temperature scale in spin dynamics. 

2. Spin susceptibility function 

The spin dynamic structure factor S(q, uj) is defined as a Fourier transformation of the 
spin-spin correlation function 

S a p(q,uj) = —jjte^ t j d 2 re-^ r < S a (r,t) • S^(0,0) > . (3.11) 

S a p(q,uj) can be directly measured in neutron-scattering experiment, and is related to the 
dynamic spin susceptibility function x Q/3 (q, cj) through the so-called fluctuation-dissipation 
theorem 

S a ?(q,uj) = - [l -e-M _1 x^(q,a;), (3.12) 

where (3 = 1/ksT, and Xa/?(9> a ') = ^ m Xa/3( < i, & + ^0 + ) is the imaginary part of the retarded 
spin susceptibility function. 

In the Matsubara representation, the transverse spin susceptibility is defined by 
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X(q, iu n ) = J dre luJnT J d 2 re~^ r < T T S + {r, t)S~(0, 0) >, (3.13) 
with uj n = 2im/ (3. In terms of a continuum version of Eq.(2.46), one may write 

< T T S + (r,r)S-(0,0) >= j f ^ e ^» r ] < T r 6+(r, r)6 T (0, 0)e$fil A "' dr > 



r)6t(0, 

(°'°) Ah. 



Qo 



x < T T 6 i (r,r)6+(0,0)e J (^> >, (3.14) 

in which Qo = (±-,±-) are the AF wave vectors. A h in Eq.(3.14) is from the nonlocal 
phase in the spin expression (2.46). In the mean-field approximation with A h replaced by 
A , Eq.(3.14) may be rewritten as 

<T T S + (r,r)S- (0,0) >= j (^e^A e 2i fo^Gl(-r,-T)Gi(v,r), (3.15) 

where the line integration on the right-hand side is chosen along a straight line connecting 
(x, y) and (0, 0) on 2D plane. The Green's function G% is defined by 

G£(r,r) = - < T r 6 CT (r,r)6+(0,0) > . (3.16) 

b a (r) may be expressed in the representation of the Landau levels as b a (r) = J2nk < r\nk > 
b° k , where n = 0, 1, 2 • • • , is the Landau level index and k is the quantum number inside 
each Landau level. Correspondingly, 

G CT (r,r)= ]T < r\nk X n'k'\0 X T r ^(r)^, fc , + (0) > 

n, k n' ,k' 

= ~ E < r l (E I"* >< 1° >< T T b° nk (T)b° nk + (0) > 

n \ k ) 

= J2K(r,0)G° b (n,T). (3.17) 

n 

In obtaining last line of Eq.(3.17), the fc-dependence of the Green's function Gb = — < 
T T b n k( T )bnk(.fy > nas been neglected. Note that k represents the center of each cyclotron 
orbital in the present symmetric gauge and ^-dependence of Gb should not be important 
due to the translational invariance of the system. IT^(r, 0) in Eq.(3.17) is given by@ 

n*(r,0) = L n ^W(r,0) (3.18) 
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with L n (t) as the Laguerre polynomials: L (t) — 1, Li(t) — 1 — t, etc., and 

W (r,0) = ^exp[-r 2 /All\. (3.19) 

Then one gets 

X (q,^ n )= E ^m(q) / — Ph {l,u')p h {m,u") } -f - 1 3.20 

Zj ^ J-oo 27T ^ n + a;' - to" 

where 



ivL,(y)L m (y)e->Mj2i& |q - Q„|), (3.21) 



167T „ 

Qo 



with Jo(x) as the Bessel function, and the spectral function pb is defined through 

C?((,^)=/^4, (3.22) 
J 27T «w n — a;' 

in which the spin index a has been omitted for simplicity since pb does not explicitly depend 
on it in the unpolarized case. n(oo) = — 1) in Eq.(3.20) is the Bose distribution 

function. Note that the spectral function pb(l,ou) ensures uo > ujq — p > 0, where ujq is 
the energy minimum of spinon spectrum and p is the chemical potentail determined by 
E/(2tt/ c 2 )- 1 f(du/2n)p(l, u)n(u) = (1 - 5) /2a 2 . 
Finally x"(q, a>) is found by 



/oo doj' 
—p b {l, u/)p b (l', oj + u') [n(u + u') - n(u')\ . (3. 
-oo ZTT 



23) 



In the next sections, we shall examine the q-dependence, ^-dependence and temperature- 
dependence of x"(q, a;) and compare them with neutron-scattering and NMR experimental 
measurements in cuprates. 



3. Basic characteristics of spin dynamics in 2D 

Spin susceptibility x"(q, oo) in Eq.(3.23) will characterize the basic features of spin dy- 
namics. The real part of spin susceptibility can be also determined from x"(q, cj) through 
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the Kramers-Kronig relation. Here we shall mainly focus on x"^ 00 ) m the vicinity of the 
AF wave vector Q and low energy uj. 

As discussed in Sec.i, doping creates a Landau-level structure in the spinon spectrum. So 
the low-lying spin fluctuations are expected to be sensitive to doping. Due to the particular 
Landau-level structure, one can distinguish two types of contributions to x"(q, cu), which 
correspond to intra-Landau-level transition and inter-Landau-level transition, respectively. 
We will show below that these two processes are related to the so-called commensurate and 
incommensurate AF spin fluctuations in the present state. 

Commensurate AF fluctuation First we consider x"(q, at small ui such that only 
the intra Landau-level transition contributes. At low temperature with spinons staying in 
the LLL, the q dependence of x"( c l; a; ) [Eq.(3.23)] is solely decided by K 00 (q): 



in terms of Eq.(3.21). Eq.(3.24) shows that x"(q, u;) will be peaked at the AF wave vector 
Qo's. The corresponding spin fluctuation is known as the commensurate AF fluctuation. 
With the increase of temperature, spinons can be thermally excited to higher Landau levels 
such that K n (q), i^ 2 2(q), etc., will contribute. For example, K n (q) has the following form: 



which is still peaked at Q with a width essentially the same as K 00 (q) in Eq.(3.24). K u 
at higher I can be checked by numerical calculation and generally a Gaussian (3.24) is well 
satisfied when / is not too big. This means that the width is not sensitive to temperature. 
Thus, when the intra-level transition dominates, x"(q, cu) generally follows a Gaussian: 



around Q , where a = l// c and A(uS) = x"(Qo,^)- The Gaussian (3.26) with a width a de- 
termines a spin-spin correlation in real space: cos(Q -r) exp[— |r| 2 /£ 2 ], where the correlation 
length 




(3.24) 




(3.25) 




(3.26) 
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is in the same order of the average hole-hole distance. Namely, the doped holes break up 
the long-range AF correlation into short-range AF fragments with a length scale ~ £. 

The energy scale of this AF fluctuation will be characterized by the u dependence of 
A(oj) = x"(Qo? w ) m Eq.(3.26). Since only the intra-level transition is involved here, the 
magnetic energy scale will be basically decided by the broadenings of the Landau levels. 
The expression of x"(Qo,u) can be found as 

X "(Qo, u) = ~ J ^ £ Pbih <S)f»(l, " + + a/) - n{u')\ . (3.28) 

Here the broadening width Y s of the spectral function p b for a given Landau level is generally 
of the order 5 J, as caused by the fluctuation ^A^. The detailed broadenings will be sensitive 
to 5A h and other factors, which is to be further discussed later. Equation (3.28) has been 
obtained from Eq.(3.23) by using 

tf«'(Qo) = (3-29) 

107T 

in terms of Eq.(3.21) and the orthogonality of the Laguerre functions 

/ dyL l (y)L l ,(y)e-y = 8 ltV . (3.30) 

J 

[The contribution to Ku>(Q ) from the other three Q 's is exponentially small and has been 
neglected.] 

The typical oj dependence of x"(Q ,o;) is shown in Fig.l at various temperatures. The 
choice of the spectral function p b and other parameters in Fig.l is to describe the underdoped 
Y Bd2Cu30§£ n which will be explained in the next section. Here we mainly focus on the 
general features shown in Fig.l. Besides an overall small energy scale decided by the Landau- 
level broadening, Fig.l shows an interesting temperature effect characterized by the Bose- 
condensation temperature T*. When T > T c *, one finds x"(Q ,cj) oc u/T at small u>. The 
slope of the linear ui dependence increase with the decrease of T. However, this increase gets 
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arrested at T = T*, and when T < T* the low-energy part of x"(Qo, becomes continuously 
suppressed instead. This feature resembles a typical "spin gap" behavior. When T < T*, 
the spinons begin to condensate into the bottom state of the lowest Landau level. The 
contribution due to the transition from such a condensate state to the rest of the quantum 
states in the broadened LLL will then emerge in x"(Qo> This process effectively will map 
out the shape of the spectral function of the same Landau level, and becomes dominant with 
the decrease of temperature. In fact, at T = when all the spinons are condensed, x"(Qcb w ) 
in Eq.(3.28) is simply reduced to 

x"(Qo, ^)t=o = (j^j Pb{l = 0, uj) t=0 , (3.31) 

with n s — 1 — 8, which is directly proportional to the spectral function pb{0, uj) of the LLL. 
If pb(0,Lj) has a small Lifshitz tail before it vanishes at the low-energy end u = 0, a gap 
feature is exhibited in x"(Q ,o;) as illustrated in Fig.l. Even in the case that p b (0,u) does 
not approach to zero very fast at uj — 0, x"(Qo, could still generically show a pseudo-gap 
trend below T < T* as described above. This spin gap behavior and the temperature scale 
T* are the unique, and important features of the present spin state. 

Another interesting property of x"(q, cu) is its behavior at uj — > limit. This behavior 
can be probed by NMR measurement.il The NMR spin-lattice relaxation rate of nuclear spin 
due to the coupling to spin fluctuation described by x"( < l) a; ) i s given byH 

1 ^ T ^ A 2( q )X"(q^) 



Tl iV q UJ 



(3.32) 

where the form factor A 2 (q) is from the hyperfine coupling between nuclear spin and the spin 
fluctuation. For 63 Cw(2) nuclear spin in the cuprates, with the applied field perpendicular 
to the Cu02 plane, the form factor A 2 (q) is found to be0@lll 

^ 2 (q) 63Cu = [A± + 25(cos q x a + cos q y a)f , (3.33) 

where the hyperfine couplings A± and B are estimated as0 A±/B ~ 0.84, B ~ 3.8 x 
10~ 4 meV [These coefficients may slightly vary among YBCO and LSCO compounds]. For 
0(2) nuclear spin, one hasffiQ 
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A 2 



(q) 



2C 2 [1 + cos(q x a)} , 



(3.34) 



with C ~ 0.875 [Ref. 0. A 2 (q)|i7 in Eq.(3.34) vanishes at q = Q . Thus a combined 
measurement of l/ 63 Ti and l/ 17 Ti can provide a q-dependent information about x"(q, u;) 
at to — > 0. 

In the present framework, x"(q, cu) have been already obtained in Eq.(3.23). So by 
substituting it into the above 1/Ti expression, it is straightforward to get 



In obtaining Eqs.(3.35) and (3.36) only the intra-level transition (i.e, I = I') is involved. 
Using the same x"^^) whose (^-dependence at q = Q is illustrated in Fig.l, one finds the 
corresponding l/ 63 TiT vs. T as shown in Fig.2. In contrast to a usual Korringa rule for a 
Fermi liquid: \jT\T ~ const., Fig.2 clearly depicts a non-Korringa behavior which falls off 
like oc 1/T at high temperature. And l/ 63 TiT is peaked at the characteristic temperature 
T*, below which x"(Q ,co>) begins to develop a "spin gap" behavior. 



On the other hand, at oxygen 17 0(2) site, the coefficient Cj has to be replaced by Gf 
in Eq.(3.35): 



The difference between l/ 63 Ti and l/ 17 Ti is due to the different weight functions [Eqs.(3.33) 
and (3.34)], which pick up the commensurate magnetic fluctuation around q = Qo at G3 Cu(2) 
site, but suppresses such a contribution at 17 0(2) site. Cf in Eq. (3.37) simply vanishes 
at £ — > oo: the long-range AF limit where the whole contribution is from q = Qo- At 
£ = 4a, for example, we find a reduction of ~ 1/93 for l/ 17 Tj as compared to l/ 63 ^. Thus 




(3.35) 



where 




(3.36) 




(3.37) 
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the commensurate AF fluctuation as described by Eq.(3.26) has negligible contribution to 
the spin relaxation for planar oxygen nuclei of the cuprates. In contrast, non-magnetic 
incoherent contributions should not be suppressed so strongly by the form factor, and would 
become dominant contribution in l/ 17 Ti as well as in the Knight shift,i which measures the 
real part of the spin susceptibility near q = 0. 

Hence, u)/uj at uo — > gives rise to a non-Korringa behavior of the NMR spin 

relaxation rate. This 1/T law of 1/T\T can be obtained analytically, if the broadening T° s 
of the spectral function pb for the LLL is small as compared to temperature. That is, when 
r° <C T <C uj c (and of course T > T*), one has 

f3 38) 

63 TiT ~~ T ' 

in terms of Eq.(3.35), where D = ttCqUs/T^S 2 . Next let us consider the case when tempera- 
ture is further increased such that more Landau levels are involved. Here we have to assume 
a general Lorentz-like broadening for each Landau level with T s = 0Auj c oc 5. Then l/ 63 Ti 
at high temperature is shown in Fig.3, where the curves become very flat and not sensitive 
to the doping concentration. This is in contrast to the low temperature regime, where a 
strong doping-dependence can emerge (in Eq.(3.38), D ~ 1/5 if one simply takes r° oc 5). 

Therefore, a full picture for the low- lying commensurate AF correlation is formed in 
the present spin-charge separation state. This picture is rather unique, and is drastically 
different from those of Fermi liquid and local antiferromagnetic descriptions. Doping effect 
plays a key role here. It decides a doping-dependent energy scale for the AF fluctuation, 
which can be much smaller than the characteristic energies in a Fermi liquid (e/) and an 
effective local antiferromagnet ([1 — 5] J). It also leads to a doping-dependent correlation 
length, comparable with the average spacing of holes. Furthermore, the Bose condensation 
of spinons determines a new characteristic temperature scale, below which the low-lying 
spin fluctuation is suppressed and non-Korringa behavior of the spin relaxation rate gets 
interrupted, similar to a spin gap effect. 

Incommensurate AF fluctuation Now let us consider higher energy regime. If u is 
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increased such that the inter Landau-level transition gets involved, the simple Gaussian- 
like q-dependence in Eq.(3.26) will be modified. A typical x"(q, oS) vs. q is shown in 
Fig. 4 at different lu's. For the purpose of illustration, we have chosen the same broadening 
T s = 0Auj c = OA5u>o for each Landau level, and T = O.Iujq, 5 = 0.15. It shows that the 
width of the Gaussian is broadened at first with the increase of u. Then the commensurate 
peak at Qo is split and two incommensurate peaks emerge at some fixed positions as the 
inter- level transition becomes dominant. In fact, this incommensurate structure is due to 
-fToi(q) in Eq.(3.23). According to Eq.(3.21), one finds 

~ t s ^ \- |q- Qol 2 / |q-QoP\ , Q qm 

^ 0l(q) = T6^^^^ eXP ( ) ' (3 ' 39) 

which leads to incommensurate peaks at a ring circling Q by a radius of \ / 2n5/a. We 
stress that the doping concentration is presumably small here so that the lattice effect is 
negligible. At a finite doping where spin-correlation length becomes comparable with the 
lattice spacing, the lattice effect is expected to become important, which could strongly 
affect the positions of the incommensurate peaks in q space. In this case, both flux and 
lattice have to be treated on the same footing like in the incommensurate-flux-phase.il 

With the further increase of u, more complicated structure will show up as Kq2, i^03, 
etc., are involved. The weight of x"(q) will be further shifted towards whole Brillouin zone, 
i.e., non-magnetic regime. At such a high energy, the present continuum approximation 
may no longer be appropriate. Nevertheless, one could still get some generic feature. By 
integrating q within the whole q space, we get 

q 



5 2 rcLu' 
8tt 7 27 



{n{u + u')-n(u')). (3.40) 



Equation (3.40) predicts that at high energy, the q-integrated Xtotaii. 00 ) w ^ be saturated at a 
constant level, modulated by a Landau-level-like oscillation. We note that the higher energy 
part is contributed by those fluctuations with less magnetic character. Thus it is less easy 
for elastic neutron-scattering to fully collect data for x'Ltaii 10 )- On the other hand, at low 
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energy, finite widths of the Landau levels will lead to a more broadened first peak in x'lotaii 00 ) 
as compared to the overall energy peak of x"(Q ,u;) where only intra-level transitions are 
involved. 

C. Comparison with experimental measurements in cuprates 

The normal-state spin dynamics represents an important characteristics for the high-T c 
cuprates. A lot of measurements have been done on these materials, especially YBCO 
and LSCO compounds. Experimental results have revealed a rather rich phenomenon, and 
also seem to suggest that material-dependent effects, like double-layer structure in YBCO, 
may play major roles in the spin dynamics. Thus, whether there is an underlying universal 
mechanism for spin dynamics in the cuprates is not transparent in terms of the experiments 
alone. In the following, we argue that the present spin-charge separation scheme can relate 
a variety of anomalous spin properties together, and provide a consistent picture for these 
materials. 

In the underdoped cuprates, like Y Ba2Cu?,0§& and insulating La1.95Bao.5CuO 4, a com- 
mensurate AF fluctuation has been verified by neutron scattering,0'00 which is peaked 
at Q in x" ■ These neutron data have been well fitted00 by a Gaussian form (3.26), and 
the width a is indeed roughly independent of temperature. The corresponding spin-spin 
correlation length follows a a/\/5 rule similar to Eq.(3.27), as shownEl in LSCO at small 
doping. Thus the momentum feature of the commensurate AF fluctuation at small doping 
is well described by the theory. In the underdoped YBa 2 Cu30e+ x , the energy scale of such 
a commensurate AF fluctuation has been systematically investigated. Emm It is found to 
be doping dependent and small as compared to the exchange energy J ~ 120 meV. As 
discussed in the last section, a small, doping-sensitive energy scale for the commensurate 
AF fluctuation is one of intrinsic features of the present scheme, and is distinguished from 
the usual Fermi liquid as well as the local spin descriptions. To our knowledge, so far there 
is no other alternative theory could be able to obtain such a small energy scale (<^C J) for the 
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commensurate AF fluctuation. Furthermore, a prominent feature has been exhibited^ at 
the low energy part of the spectroscopy, where the weight of x" is continuously suppressed 
below some characteristic temperature, resembling an opening of a spin gap. In the present 
theory, even though there is no real gap in the spin fluctuation spectrum, a characteristic 
temperature T* is found to be naturally associated with a spin-gap-like phenomenon. In 
Fig. 1, a typical energy and temperature dependence of x"(Qo; w ) has been shown, where 
the spectral function is chosen to describe YBd2CuzO§& (see below). And the overall energy 
and temperature features in Fig.l, particularly the "spin gap" behavior, are in good agree- 
ment with those found in the underdoped yi^C^Oe+zBS In the underdoped cuprates, 
the NMR spin relaxation rate of planar copper nucleiS manifests a non-Korringa behavior, 
which is suppressed when the spin gap feature shows up in the neutron scattering at low 
temperature. In contrast, a conventional Korringa temperature behavioral is found for pla- 
nar oxygen nuclei. The theoretical 1/ 63 T\T has been presented in Fig. 2, which is calculated 
by using the same spectral function as used in Fig. 1 for YBa>2Cu 3 6+x . Fig. 2 shows a 
non-Korringa behavior at T > T* as well as a "spin gap" feature below T*. All of them 
are also qualitatively consistent with the experimental measurements [the contribution of 
the present commensurate AF fluctuation to l/ 17 TiT is dramatically suppressed (~ 1/93 
at £ ~ 4a for Y Ba2Cu 3 OQ,^) so that the non-Korringa signal does not leak to the oxygen 
sites]. Therefore, the q-, to- and T- dependence of the present low-lying magnetic fluctua- 
tion at small doping are all in agreement with the main experimental features found in the 
underdoped cuprates. 

Techniquely, let us discuss the spectral function used in the theoretical calculation. Ac- 
cording to Eq.(3.31), one may use the experimental measurement of x"(Qo,i^) at low tem- 
perature to determine p&(0,co>), instead of a first-principal calculation which would involve 
much complicate factors here like inter-layer coupling. As pb(0,u) satisfies the normalized 
condition, if its broadening at low temperature is quite large as compared to the temperature 
scale which we are interested, the temperature dependence of p b (0,u) may not be impor- 
tant. Then x"(Qo,co>) in a whole temperature range can be determined. In Fig.l, we have 
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chosen p&(0, uj) such that to give a right energy scale of x' / (Qo) w ) f° r Y Ba 2 Cu^O%,Q at low- 
temperature [Fig. 7 in Ref. pOj] . The contribution from the higher Landau levels are neglected 
because the temperature range considered is comparatively smaller than the broadening of 
the LLL. Besides the shape of pb(0,u), an overall strength of x"(Qo,co>) has also been ad- 
justed, in order to compare with experiment, by reducing the spinon concentration from 
n s = 1 — 5 to n* < n s . This is because in the present approximation the spinons are treated 
as ideal bosons instead of hard-core bosons. At half-filling, this approximation could lead to 
a too big magnetization. In the present doped case, it would also give rise to a too strong 
magnetic fluctuation. The hard-core effect as well as the incoherent band in the spinon 
spectrum should reduce the effective number of spinons contributing to the AF correlation. 
We find a reduction of n*/n s ~ 1/3 at 5 = 0.10 giving a T* ~ 175K, in consistence with 
the corresponding experimental characteristic temperature (~ 160 K), and at the same 
time, leading to a spin relaxation rate whose magnitude quantitatively agreeing with the 
NMR measurement .Si 

To end discussion of the underdoped materials, we give several remarks below. Some 
authors&H have attributed the spin gap phenomenon in the underdoped YBCO to their 
peculiar bi- layer structure. In contrast, one may have noticed that the theory here is purely 
two dimensional. Nevertheless, it has been noted that the spinon spectral function pb(0,uj) 
for the LLL has been determined directly from the experiments, in which the inter-layer 
coupling could have been already included and may play a key role in the broadening. 
Recall that in the present scheme, there is no real gap opened in the spin spectrum, and 
the "gap" feature is very sensitive to the detailed broadening of the LLL. Theoretically, 
for a pure 2D system one would expect the broadening of the LLL to be much narrower 
at low temperature. The reason is that a spinon in the LLL could not be scattered into 
a lower energy state while emits a phonon-like excitation of SA , because it is already at 
the energy bottom. It cannot jump up to higher Landau-levels either, due to the energy 
conservation. Such a feature has been indeed found in a similar problem.0 If this were true, 
the energy scale in the underdoped material might have been much sharper than observed. 
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But in the insulating cuprates like Lai^Bao^CuO^ [Ref. [17]], holons should be localized and 
randomly distributed such that 5A h describes a random flux with energy scale ~ 0, instead 
of a phonon-like dynamic mode. Then the degeneracy of the LLL can be lifted under the 
static random vector potential, and consequently a much broadened LLL could appear. 
And mixing with other levels is also expected here. For the underdoped YBCO, however, 
one is in the metallic phase, and holes must be mobile which would lead to a well-defined 
dynamics for 5A h . Nonetheless, one may still expect a strong localization tendency of holes 
at low temperature, as suggested by the transport measurement s J3 which could in turn lead 
to a more broadened LLL than in the optimally-doped regime. Furthermore, the bi-layer 
coupling in YBCO can split a sharp LLL into two peaks (symmetric and antisymmetric 
states), but for a larger broadened LLL, the bi-layer coupling may well result in a single, 
much broadened peak like the one shownS in Y Ba^Cu^Og^. Finally, we point out that in 
these underdoped materials, incommensurate structure has not been observed yet at high 
energy. But the width a in Eq.(3.26) has been found to be increased with the energy, which 
is an indication of the involvement of inter-level transition in the present theory. 

Next we consider optimally doped cuprates, like YBa2Cu 3 0j^ y and metallic 
La2- x Sr x CuO\. A striking common feature in these materials is the lack of a commen- 
surate AF fluctuation at low energy in neutron measurements. The spin-polarized neutron- 
scattering measurement^ of Y Ba 2 Cu- i O-j^y has only revealed a sharp high energy AF peak 
around 41 meV. In the metallic LSCO (x = 0.075, 0.14 and 0.15), an incommensurate 
structure has been foundEl down to an energy scale ~ 1 meV, and no commensurate AF 
correlation is observed within the experimental resolution. However, the absence of a low- 
energy AF fluctuation or the presence of a low-energy incommensurate fluctuation would be 
both in sharp conflict with the NMR measurements^ which probes u ~ regime and implies 
a strong commensurate AF correlations in these materials. For example, by extrapolating 
the spin susceptibility with the incommensurate structure observed by neutron-scattering 
down to the NMR frequency (~ 10~ 4 meV), it has been found!!! that a large magnetic 
contribution could leak to 17 0(2) sites, in contrast to the measured spin relaxation rate 
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l/ 17 Ti which appearsB to be completely dominated by a nonmagnetic (~ q independent) 
contribution. As a matter of fact, the spin relaxation rates in LSCO have shown canonical 
behaviors, which basically are the same as those found in YBCO compounds. As elaborated 
in Sec. B3, only a commensurate contribution can lead to a non-Korringa behavior of 1/Ti at 
63 Cw(2) sites, and at the same time, be strongly cancelled out at 17 0(2) sites. This fact has 
been the basis for the so-called nearly antiferromagnetic-Fermi-liquid theoryffl Therefore, 
in order to reconcile neutron and NMR experiments, one is led to the conclusion that a 
commensurate AF fluctuation should re-emerge in the optimally-doped cuprates within an 
energy scale beyond the experimental resolution. 

However, it is hard to perceive such a characteristic scale in a conventional theory. So 
far, there have been a number of theoretical conjectures for the mechanism of incommensu- 
rate magnetic fluctuations in the metallic LSCO. In the Fermi-liquid-like framework, the 
incommensurability is directly connectedS0 with the Fermi-surface shape. But it lacks a 
small energy scale, within which a commensurate structure could be recovered. The spiral 
state of Shraiman and SiggiaS based on the t — J model also provides an incommensurate 
structure. The original spiral state is a long-ranged state, and some short-range versions 
have been proposed.BS But a low-energy commensurability is still hard to be comprehended 
here. Some non-intrinsic mechanism for the incommensurability due to the inhomogeneity 
in the LSCO system has been also proposedi! in order to reconcile the neutron and NMR 
data, where the commensurate AF fluctuation is assumed intrinsic in the metallic LSCO. 
Nevertheless, the details still need to be carried out in order to make comparison with the 
experiments. 

The present spin-charge separation scheme is unique to have a small characteristic energy 
scale of the commensurate AF fluctuation, as already discussed in the underdoped case. As 
pointed out there, an even smaller energy scale could be present at optimal regime as the 
LLL broadening becomes sharper when holons are much mobile in this larger doping regime. 
Due to the narrowness of the energy range, the amplitude of x" will be also small due to the 
cancellation of the Bose functions inside Eq.(3.23) [one has x" ~ > at uo/T — > 0]. Therefore, 
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this commensurate AF correlation may well be beyond the experimental resolution to be 
directly observed by neutron-scattering. In YBa 2 Cu307_ y , the bi-layer coupling could also 
split the sharp LLL into two peaks as mentioned before. It then explains the sharp 41 meV 
AF peak found by neutron-scattering0 as a result of the spinon transition between these two 
split peaks. On the other hand, such an AF fluctuation can give rise to a non-Korringa law of 
the NMR spin relaxation rate [Eq.(3.38)], which is consistent with the NMR measurements 
in YBCO and LSCO. Imai et a/S have measured l/ 63 Ti in La2- x Sr x Cu?Q \ system up 
to 900K. At high temperature, all data (x = — 0.15) seem to converge and saturate to 
the same temperature-independent value 2700 ± 150sec _1 . The theoretical results shown in 
Fig. 3 agree well with this tendency. And if we choose, say, u = lOOO-ft', we find that the 
saturation value of l/ 63 Ti is about twice larger than the measured one, which is a reasonable 
value, considering no adjustment of the spectral function has been made to fit the data. 

In our theory, the incommensurate fluctuation in the metallic LSCO will be attributed 
to the dominance of the inter-Landau-level transition in the experimental energy-transfer 
regime, as discussed in Sec. B3. To be consistent with neutron-scattering measurements, the 
LLL broadening has to be very sharp as explained above, while the second level broadening 
is relatively larger, which should be centered around u c ~ 10 meV in x = 0.15 case. When 
the temperature is increased such that ksT ~ u c , a large amount of spinons are expected 
to be thermally excited to the second Landau level. Then in the experimental-observable 
energy regime, the intra-level transition could emerge again due to a larger broadening in 
the second Landau level. Correspondingly, the incommensurate structure should be replaced 
once again by the commensurate peak around Qo- This has been indeed observed in the 
neutron-scattering measurement Jll where a broad commensurate peak is found to re-emerge 
around T ~ 100AT. 

The characteristic energy scale uj c of YBCO seems to be several times larger than that of 
LSCO for some unknown reason. This decides a main distinction of spin dynamics between 
YBCO and LSCO, in terms of the present theory. In other words, one should expect that at 
a sufficient high energy, incommensurate structures could also show up in YBCO system. Of 
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course, many factors, especially the inter-layer coupling, may complicate the details. Due to 
these incommensurate contributions at high energy, the q-integrated susceptibility function 
Xtotai [Eq.(3.40)] will stretch up to an overall energy scale ~ J s with a roughly-constant 
amplitude (modulated by the oscillation), Such a unique Xtotai behavior is a high energy 
prediction of the present theory for the cuprates. 

In the optimally-doped LSCO and YBCO, the spin gap feature is absent in the com- 
mensurate spectrum due to the fact that the energy scale is too small. But according to 
the theory, a pseudo-gap should be also present in the LSCO compounds below T*, which 
involves the inter-level transition instead of intra-level transition in the commensurate fluctu- 
ation case. This pseudo-gap phenomenon has been clearly shown by the neutron scattering.!^! 
In these optimally-doped cases, the spin characteristic temperature T* should become very 
close to the superconducting transition temperature T c , to be consistent with neutron scat- 
tering and NMR measurements. In the present theory, superconducting transition will occur 
when spinons and holons are both condensed. Usually one finds T* > T c at small doping. 
But when the spinons are condensed, the frustration effect on holons from the spin part 
(see next section) will be reduced too, which in turn is in favor of the Bose condensation of 
holons. In other words, T c and T* may correlate with each other. An optimal regime in the 
present theory may be properly defined as when T* and T c coincides. A further discussion 
of the superconducting transition will be presented elsewhere. 

IV. TRANSPORT PROPERTIES 

Due to spin-charge separation, holon degree of freedom is to be solely responsible for the 
electron transport phenomenon in this system. As noted before, the nonlocal phase in 
the holon Hamiltonian Eq.(1.3) will provide an unconventional scattering mechanism in two 
dimensions, as in contrast to the ID case where it vanishes and holons simply behave like free 
particles. In the following, we will explore 2D transport properties under such an effective 
Hamiltonian H h . Similar to the magnetic properties discussed in Sec. Ill, transport will be 
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an another crucial test for the experimental relevance of the present spin-charge separation 
theory. 

For simplicity, we may consider Hh in the continuum limit at small doping. The validity 
of the continuum approximation will be discussed later. The corresponding continuum 
version of Eq.(1.3) has the following foruS 



where the holon mass m t = (V^tha 2 )- 1 , and A s is defined by A s = A+ + A!, with 



Here p s a {r) = b^(r)b a (r) is the local spinon density with spin index a. Similar to the lattice 
version [Eq.(1.5)], equation (4.2) describes fictitious 7r-flux quanta bound to spinons, and 
the sign a = ±1 in front of the integration means that spinons with different spins carry flux 
tubes in opposite directions (see Fig. 5). In an unpolarized system, one has < pi >=< pj > 
so that on average < A s >= 0. In other words, the scattering source as contributed by 
spinons may be regarded as a sort of fluctuating gauge field, and no T- and P-symmetry 
violations occur here. 

Gauge-fluctuation-related scattering mechanisms have been intensively studiedly 
within the framework of the uniform RVB state. So we need to distinguish the present 
gauge field A s and those studied in the literature. Since the fictitious flux quanta are bound 
to spinons in the present case, let us first consider the local density distribution of spinons. 
At small doping, spinons are forming fluctuating AF domains, as discussed in Sec. Ill, within 
the spin correlation length £ > a. For the purpose of illustration, an extreme case is shown 
in Fig. 6, where £ ~ a and only a pair of f and J. spinons are bound together like a valence 
bond. + and — signs in Fig. 6 represent a spin configuration, or equivalently, directions of 
the 7r-flux quanta attached to spinons. If a holon moves about a closed path C as shown 
in Fig. 6, the accumulated flux will obviously have a large cancellation within the loop, with 
only the contribution from those near the path C whose pairs are cut by the loop. Then the 




(4.1) 




(4.2) 



51 



mean-square accumulated flux will follow a perimeter law: oc Lc {Lc is the perimeter of the 
loop C). It is easy to observe that this perimeter law is generally true as long as there exists 
a short-range AF correlation, which is present at small doping and within a temperature 
range ksT < J. This is in contrast to the gauge fluctuation in the uniform RVB state where 
fluctuating flux satisfiesilS an area law instead of a perimeter law. And the present flux 
problem is more similar to a random flux problem in the so-called "Meissner phase" ,@ where 
a perimeter law is present. But in the conventional case, the strength of the fluctuating flux 
is presumably weak so that the accumulated flux enclosed in a loop will become negligible 
when the loop is small enough. In this case only the long-distance behavior matters, and 
a perturbative approach may be applicable which is similar to a usual phonon-scattering 
problem. 

However, in the present case, each flux tube attached to a spinon is quantized at tt. It 
means that a slight deformation of the path C, with one spinon enclosed or excluded, could 
lead to an additional Berry-phase ±7r, or a sign change of the wavefunction. A strong phase 
interference is therefore expected at a short- distanc^ due to the high density of spinons at 
small doping. Such a short-distance effect can drastically change the nature of scattering 
mechanism which in the usual case would only involve long-wavelength processes. This 
is going to be a key distinction between the present theory and the usual gauge theories. 
We will see that such a short-range interference can lead to an exotic "localization" effect 
of holons, and result in a set of very interesting transport anomalous in longitudinal and 
transverse channels. 

We begin by considering a single holon problem. For an one-body problem, a Feynman 
path-integral formulation will become very useful. The transition amplitude for a holon to 
travel from a to b in space under the vector potential A s is given by0 




(4.3) 



with the action 




(4.4) 
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By introducing the following identities 

W) 



^ 2 - e i ^ M «) = tf(r-P 9 ), (4.5b) 



K(b, a) may be rewritten as 

K(b,a) = J° 1 Vv{t)Vr p {t)Vv q {t) J VX (t) J V(3 (t) e iS ^ [b ' a] , (4.6) 



in which 



S xp [b, a]=f t b dt [^r 2 + ^ + ^r 2 ] - / dr p • A ;(r p ) - | dr q • Ai(r„ 



+ / dt[A(t) • (r - r p ) + (3(t) • (r - r,)], (4.7) 

where + m p + m q = m t . Extra degrees of freedom are thus introduced in Eq.(4.6). 
Without the last term involving the Lagrangian multipliers A and (3, S\p[b, a] in Eq.(4.7) 
would simply describe three independent particles: a free holon with mass ra^ known as h 
species, p and q species which interact with the vector potentials A+ and At, respectively. 
The fields A and (3 play the role to recombine these three species together as a real holon, 
and thus effectively eliminate the additional degrees of freedom in the end. 

So far no approximation has been made. By introducing p and q degrees of freedom, 
the effects of A^_ and A s _ are separated. In terms of Eq.(4.2), A* may be rewritten as 
A* = a X s + 5A s a , where 

A s = ^|(zxr), (4.8) 

with B s = (1 — 8)tt/2, and 5 A* has the same form as Eq.(4.2) but with p s a replaced by 
Sp s a = p s a — < p s a >. Similar procedure has been used in dealing with the topological phase 
A h in spinon system (cf. Sec. Ill), and is familiar in an anyon problem.ili! Here oA s 
describes a mean-field effect with the flux-quanta smeared out uniformly in space. And 5A s a 
represents the fluctuation of the flux-quanta due to the density fluctuation of spinons with 
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spin index a. In the following, we shall show that, to a leading order of approximation, 
the effect of 5A s a in the action (4.7) could be effectively represented by a relaxation of the 
binding constraint enforced by A and f3 within a scale of the cyclotron length d c of p and 
q species (under the fictitious field B s ). Under such an approximation a renormalized holon 
will be a composite particle of a bare holon and the p and q species which are in the cyclotron 
orbitals. Below we elaborate this approximation. 

Due to the symmetry, one only needs to focus on the p species. Let us consider an 
arbitrary closed path C on a 2D plane (Fig. 7). The contribution of A+ to the transition 
amplitude in Eq. (4.6) for such a path will be a gauge-invariant phase: 

f c dr p - A; = f c dr P - A s + ^dr p ■ 5A S + . (4.9) 

Eq. (4.9) may be further rewritten as 



dr v ■ A s 

c 



I dr p -A s , (4.10) 

J Cr> 



in which a path C p is introduced such that 

= ttAJVJ. (4.11) 

Here the path C p as shown in Fig. 7 is a deformation of the path C to account for the 
fluctuation of 5Af in terms of the mean- field A s . In second line of Eq.(4.11), AiVj represents 
the total fluctuating j spinon number (with regard to the average one) enclosed by the path 
C, and such a term describes the total fluctuating fluxes bound to f spinons as deduced 
from the right-hand-side integral in the first line. So the flux enclosed between the paths C 
and C p under a uniform field B s will represent the fluctuation of the flux quanta attached to 
| spinons inside the path C. In the following we shall point out that 7rA7Vj should satisfy 
a perimeter law discussed at the beginning of this section. 

In the spin background, there is no density fluctuation due to the no-double-occupancy 
constraint, and locally an increase of j spinons is always compensated by a decrease of j 
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spinons. Thus the distribution of | spinons will actually reflect that of all spins in the hole- 
absent region. With the presence of a short-range AF correlation (the correlation length 
£ > a) , excess and deficit of j spinons are neighboring to each other (a £ ~ a case is shown 
in Fig. 6). In such a "Meissner phase", the net contribution to AiVj mainly comes from j 
spinons close to the path C, and AiVj should satisfy a perimeter law instead of an area 
law as explained before. Correspondingly, the path C p should be always near the path C 
to account for AiVj in terms of Eq. (4.11), as shown in Fig. 7. An average separation d p of 
C p from C may be estimated as follows: nd^ x B s ~ one flux quanta = 7r, which leads to 
dp ~ 1/ \fB~ s — d c — the cyclotron length under B s . 

So after one replaces A+ by A s in Eq.(4.7), the fluctuation effect can be taken into 
account by deforming the path of p species from C to C p . The deviation of C p from C 
depends on the detailed fluctuation of A^_. But if temperature is sufficiently higher than 
the spinon characteristic energy scale (the broadening of the Landau level in the spinon 
spectrum), the path C p shown in Fig. 7 may be reasonably regarded as a random one, with 
an average separation d c from the path C. This effectively means that the binding constraints 
of h, p and q implemented by the Lagrangian fields A and f3 have been relaxed within a 
scale ~ d c . Here the kinetic energy cost of the deformation path C p is neglected due to the 
smallness of d c . In the weak gauge-fluctuation case, d c could be too large for the present 
approximation to be valid. 

Thus, the motion of a holon under the influence of the gauge field A s can be effectively 
described as a bare holon bound to a pair of auxiliary species, p and q, which are undergoing 
cyclotron motions in opposite directions. These p and q species reflect a "localized" effect 
caused by the phase interference at short-distances. As p and q are confined in the Landau 
levels, a large degeneracy is involved here. In real space, such a renormalized holon would 
look like a polaron, and behave like a diffusive particle in the absence of external fields. This 
peculiar structure will decide a unique transport phenomenon. 

The generalization of the above scheme to a many-body case is straightforward. Here 
one has to be cautious about the statistics of each species in the many-body case. A holon as 
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a composite particles of h, p, and q species has to satisfy the ( hard-core) bosonic statistics. 
Thus a symmetric choice would be that the h species corresponds to a boson while the p 
and q species are spinless fermions (so that the hard-core condition can be automatically 
realized). Then the effective many-body Lagrangian as a generalization of the one-body 
approximation discussed above can be written down in the functional-integral formalism as 

C h = J d 2 r [h + d T h + p + d T p + q + d T q + X(p + p - h + h) + p(q + q - h + h)) 

+ / dh W) HV + aeXtf h{r) + p+(r) HV ~ AS)2 p(r) + g+(r) + ^ - g(r) } , (4.12) 
j [ 2m/i 2m p 2m g J 

in which a e3rf is an external vector potential. And the fields h(r), p(r) and g(r) are in the 
coherent representation. The Lagrangian multipliers A(r) and f3(r) enforce the following 
binding constraint 

h + (r)h(r) = p + (r)p(r) = q + (r)q{r), (4.13) 

in the length scale larger than d c . An ultraviolet momentum cut-off A = l/d c is then implied 
in A and (3 fields in the Lagrangian L^. The masses m^,m p and m q are determined up to 
that their total summation is equal to m t . The value of each individual mass has to be 
decided beyond the present approximation, and this uncertainty will not affect the general 
conclusions we shall draw from the Lagrangian C^. 

Therefore, we obtain an effective long- wavelength Lagrangian (4.12) for the holons after 
the short-distance phase-interference is carefully treated. This is a rather unusual Lagrangian 
because there are three auxiliary species involved. But it is not derived in its first time here. 
A similar one has already been found@ in the Scheme One flux-binding state (see Sec. II) 
through different method. Even though the approximations involved in these two states are 
different, their origins are quite similar and it is not of surprise finding their normal states 
to be so close. In fact, the structures of two effective Lagrangians for holons are identical, 



except that B s is twice larger in Ref. 173 and the h field is a fermionic one there. A larger B s 



in Ref. [74] is due to the 7r-flux phase 0? has been incorporated into B s at continuum limit. 



These differences do not change the canonical behaviors described below. 
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The effective Lagrangian Ch determines an anomalous transport phenomenon!! which 
amazingly matches the essential characteristics of the optimally doped cuprates. We shall 
outline the main results in the following. For detailed discussions, one is referred to Ref. 



71 . Ch can be treated by the standard gauge-theory approach. Here it is even simpler for 
lacking the transverse fields. It is easy to show that the longitudinal fields A and (3 will 
enforce the following current constraint among h, p, and q species: 

4 = 4 = 4, ( 4 - 14 ) 

where the superscript I implies the longitudinal channel. This constraint is consistent with 
the density constraint [Eq.(4.13)]. It is important to note that there is no similar current 
constraint in the transverse channel due to the absence of the transverse gauge fields. The 
total response to an external electromagnetic field will be related to each species through 
Eq.(4.14). And different combination rules will thus be found in the longitudinal and trans- 
verse transport channels, which will lead to distinctive Hall angle behavior. 

The scattering rates of h, p and q are decided by their coupling with the longitudinal 
gauge fields A and (3 in Eq.(4.12), whose dynamics are in turn determined by coupling to h, p, 
and q species. A self-consistent treatment is required here. Due to the peculiar feature that p 
and q stay in the Landau level, one finds that the scattering rate for p and q goes like linear in 
temperature, i.e., ^- = 2KksT, where k ~ 0(1) is independent of the coupling constant (the 
masses) and only has a weak doping-dependence (and ~ oc u is also found when hu > ksT). 
And for the h species, the scattering rate has a T 2 behavior: ^- oc ^p~- m terms of the 
constraint Eq.(4.14), the total longitudinal resistivity can be determineclB and at k B T <C t h 
it is dominant by l/r s so that p oc T. The linear-T resistivity in the cuprates is indeed 
found! to be related to a linear-T relaxation rate. Particularly, the coefficient is roughly 
around a numerical factor of 2 for all the optimal materials.i This is very interesting feature 
and is quantitatively consistent with the present theory, where the coupling-independent 
coefficient 2k, ~ 2 is determined by the unique structure in the scheme. The transverse 
resistivity p xy can be also obtained, and as noted before the distinctive combination rules in 
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the longitudinal and transverse channels will lead to new consequences. Namely, the Hall 
angle as defined by cot = p/p yx is found to be related to the second scattering rate 
l/rh- cot oc 1/Th. oc T 2 . Consequently the Hall coefficient Rh follows a 1/T behavior. 
The involvement of a second scattering rate oc T 2 in the Hall angle for high-T c cuprates was 
first pointed out by Princeton group [Ref. [||| based on the analysis of the experimental 
data. The effective Lagrangian in Eq.(4.12) provides a microscopic theory for it in the first 
time. A magneto-resistance with Ap/p ~ T~ 4 dependence has been also predicted^ in the 
longitudinal channel, which has been recently observed^ in YBCO (We note that the overall 
sign is uncertain in the theory while experimentally it is found to be positive). Due to the 
"localized" effect of holon, a strongly doping-dependent thermopower has been obtained in 
the present framework, which also agrees well with the experimental measurements in the 
high-T c cuprates. 

The transport properties determined by the effective Lagrangian Ch may be called as 
the canonical ones, which well account for the optimally doped cuprates. In the following 
we briefly discuss the condition for deviation from such a canonical case in the present 
theory. The key assumption involved in deriving Eq.(4.12) from is the randomness of 
the flux quanta movement near the path C shown in Fig. 7, which leads to uncorrelated 
segments along the path C v so that p and q can be effectively treated as detached from h 
at a scale of d c . In this case, the detailed spin dynamics becomes irrelevant. In the optimal 
doped cuprates, the spinon LLL broadening is very sharp (cf. Sec. Ill) so that the above 
condition may be always satisfied at the normal-state temperature. For the underdoped 
YBCO, however, the LLL broadening is quite large, and when the thermal energy is less 
than it the detailed spin dynamics could get involved in C p , especially when T ~ T*, and 
spinons begin to condensate into the energy bottom. In this case, C p may not be treated 
as a random one with regard to C, and the present approximation could break down. As 
a consequence, a deviation from the canonical transport behaviors is expected below some 
temperature scale, which should be correlated with T*. Experimentally, such a correlation 
between the transport and the "spin gap" is indeed foundl^Jli^l in the underdoped YBCO 
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system. Furthermore, at sufficiently small doping, the transport properties described by 



Eq.(4.12) itself can also deviate from the canonical behavior as discussed in Ref. WA 



Finally, we would like to make a comment on the continuum approximation of Hh- At 
small doping, the continuum approximation can be well justified for the spinon Hamiltonian 
H s because in it is vanishingly small. In the present holon case, however, the strong 
fluctuation effect at short- distance is very important, even though < A\- >= 0. Thus one 
would expect the lattice effect to be involved even at small doping. So far we have not been 
able to include such a lattice effect. Nevertheless, as the phase interference at short- distances 
is not a coherent effect, the role of the lattice effect may not be really crucial. Furthermore, 
p and q are just auxiliary particles to take care of the phases interference, which can be 
introduced in the continuum space even when lattice is included. A further study is still 
needed. 



V. CONCLUSION 

In this paper, we have obtained the spin-charge separation scheme based on a saddle- 
point state of the t — J model. In this saddle-point state, we find a deconfinement of spin and 
charge degree of freedom at finite doping in 2D case, where the transverse gauge fluctuation 
as the confinement force is gapped. Such a gap disappears at half-filling, where spinons 
are presumably confined to form spin-1 excitations, and a long-range AF order has been 
recovered. This saddle-point state has been also shown to reproduce the known asymptotic 
spin-spin correlation in one dimension at both half-filling and finite-doping. Thus, in some 
important limits where the behaviors of the t — J model are known, the present state has 
produced the right results. These constitute important check for a strongly-correlated model 
where the conventional approximation breaks down. 

The most interesting properties have been found at finite doping in 2D case, when spin 
and charge become deconfined. In this regime, the saddle-point state becomes meaningful 
due to the suppression of the gauge fluctuation, and spinons and holons can be appropriately 
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treated as separated systems in terms of conventional approaches. We have shown that in 
this saddle-point state there still exist some exotic residual couplings between spin and 
charge degrees of freedom, in spite of the spin-charge deconfmement. These couplings are 
nonlocal in the sense that spinons can feel the existence of holons nonlocally by seeing some 
fictitious flux-quanta bound to the latter, and vice versa. Different from a usual electron- 
phonon system where the coupling is a single interactive term, spinons and holons here 
are scattered by each other in distinctive forms. These scattering forces lead to anomalous 
spin dynamics and transport properties in the present system. For example, a sharp AF 
peak centered at Q is exhibited in the imaginary dynamic spin susceptibility with an 
energy scale much smaller than the exchange energy J. The width of such an AF peak 
in q-space is determined by doping concentration in a form oc \^6, which is temperature- 
independent but increases with the energy transfer. The NMR spin relaxation rate of nuclear 
spin due to coupling with such a magnetic fluctuation shows a non-Korringa behavior for 
planar Cu nuclei and a strong suppression for planar O nuclei. Furthermore, we have 
found a characteristic temperature T* below which all these AF anomalies get suppressed, 
in resemblance to a "spin gap" system if the superconducting temperature T c < T* and 
T c < T < T*. Incommensurate AF fluctuation has also been found in this system at a 
higher energy scale. For the charge degree of freedom, we have demonstrated that the 
scattering from the spinon background leads to a strong phase interference at short-distance 
for holons, and an effective long-wavelength Lagrangian is derived. Such a Lagrangian has 
been found to give the following canonical transport phenomena: resistivity p ~ T with 
H/t ~ 2k B T; the Hall angle cot Oh oc l/r h oc T 2 ; the magneto-resistance Ap/p oc T -4 ; and 
a strong doping-dependence of the thermopower. 

The magnetic and transport properties of the present spin-charge separation state share 
amazing similarities with those found in the high-T c cuprates, as discussed in the context 
of the paper. Based on the theory, a consistent picture can be conjecturally formed for the 
normal-state of the cuprates. The optimally-doped materials can be defined as T c ~ T*, 
where the commensurate AF fluctuation energy scale becomes very small, maybe unde- 
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tectable in terms of the present neutron-scattering resolution. Such a small energy scale is 
a unique feature of the present state, and as discussed in the context, NMR and neutron- 
scattering data in both metallic LSCO and YB 2 Cu 3 0j can be reconciled here. And in the 
transport channel, the canonical behaviors are exhibited in agreement with the transport 
measurements. A larger commensurate AF fluctuation energy scale can be realized at a 
smaller doping regime, where holons tend to be localized, and the interlayer coupling in 
YBCO could further help the broadening of the energy scale. In this regime, a commen- 
surate AF fluctuation may become observable in neutron-scattering like in the underdoped 
YBCO, and pseudo spin gap behaviors will also show up in both neutron-scattering and 
NMR spin relaxation rate below T* (> T c ). In this case, a deviation from the canonical 
behaviors are expected in the transport channel below a temperature scale around T*. We 
note that in the present theory, there is no direct experimental input except that the under- 
lying t — J model is widely perceived as a simplified description of the copper-oxide layer in 
the cuprates. This becomes a compelling fact as so many important experimental features 
are naturally exhibited in the present state. We have pointed out in the introduction that 
the spin-charge separation is a key for all these delicate magnetic and transport anomalies 
to appear. 

Two other important issues remain to be clarified in the present spin-charge separation 
scheme. One is about the single-electron properties, particularly the location of Fermi 
surface, and the other is about the superconducting condensation in this framework. Since 
an electron is described as composed of two bosonic holon and spinon, the accompanied 
nonlocal phase fields in Eq.(1.6) will play a central role to restore the fermionic properties 
of the electron. As outlined in the introduction, the nonlocal phase fields will be responsible 
for an electronic Fermi surface satisfying the Luttinger volume theorem as well as a finite 
pairing order parameter < c^c^ > (and its symmetry) when both holons and spinons are 
Bose condensed. We shall discuss these important problems in the follow-up papers. 
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CAPTIONS 



Fig. 1. x"(Q ,cg>) vs. to. A pseudo spin-gap behavior is exhibited at T < T*, where the 
low-energy part (~ 5meV) is continuously suppressed with the decrease of the temperature. 
The spectral function in x"(Q ,co>) is chosen to describe Y Ba2Cu%0§.<a (see context). 

Fig. 2. The NMR spin-lattice relaxation rate for the planar Cu nuclei calculated in 
terms of the same spectral function used in Fig. 1. 

Fig. 3. High-temperature behaviors of the Cu nuclear spin relaxation rate at various 
doping concentration. Notice that the spin relaxation rates are leveled off at high tempera- 
ture and essentially indistinguishable at small doping (5 < 0.2). 

Fig. 4. Typical u ) vs - <1 a ^ various u's. With the increase of a>, the Gaussian form 
first gets broadened (dotted curves with uj = 0.10 and 0.15 [u> as the unit]), and then an 
incommensurate split occurs when the inter-level transition becomes dominant (solid curves 
with uo = 0.18, 0.20 and 0.30). 

Fig. 5. Fictitious flux-tubes attached to the spinons (solid circles) which are seen by the 
holon (open circle) on a 2D plane. 

Fig. 6. An extreme case of short-range AF correlation. The + and — signs represent 
a spin configuration, or equivalently, an array of fictitious flux-tubes bound to the spinons. 
The solid closed path cuts through those dimer pairs of spins and the total flux enclosed 
satisfies a perimeter law (see Ref. |98|) . 

Fig. 7. An arbitrary closed path C of a holon (solid one) and a corresponding path C p 
for p species (dashed one). 
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